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PREFACE 


This volume is the second in a series of reports being prepared by 
Professors Kane and Likins on various aspects of spacecraft attitude 
dynamics. The first volume, "Kinematics of Rigid Bodies in Spaceflight," 
was presented as Technical Report No. 204 of the Stanford University 
Department of Applied Mechanics, published in May of 1971. Subsequent 
volumes will deal with non-gravitational forces and moments on space- 
craft, with the attitude dynamics of rigid spacecraft, and with the 
attitude dynamics of nonrigid spacecraft. In time, the series will be 
assembled as a book intended both for classroom use and as a reference 
for practicing engineers. 

Most of the development herein is the unsponsored enterprise of the 
authors, forming the basis for the spacecraft attitude dynamics studies 
of Contract NAS8-28358. 

The present volume, "Gravitational Forces and Moments on Spacecraft," 
furnishes a unique perspective on a subject that was developed in its 
classical form long before the era of modern space exploration. The 
development provided here differs from the traditional presentation of 
the subject by focusing attention on applications to problems of 
attitude dynamics of the spacecraft. In the opinion of the authors, 
understanding of gravitational forces and moments required for the 
solution of problems of attitude dynamics can be communicated most 
directly in terms of arguments based on Newton's law of gravitation, 
presented by employing the methods of Newtonian (vectorial) mechanics, 
with minimal recourse to the classical concepts of potential theory. 
Accordingly, in this report, the necessary ideas are developed and 
relationships are established to permit the representation of gravitational 
forces and moments exerted on bodies in space by other bodies, both in 
terms involving the mass distribution properties of the bodies, and in 
terms of vector operations on scalar functions classically described as 
gravitational potential functions. 
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2.1 


2*1 Gravitational interactioB of two particles 

A particle F of mass m experiences in the presence of a particle F of 
mass m a force f acting along the line joining P and directed from P 

toward P, and having a magnitude proportional to the product of m and m * 
and Inversely proportional to the square of the distance between P and P* 
Hence, if £ is the position vector of P relative to P (see Fig. 2.1.1), 
the force F can be expressed as 

^ » ~Gnmi£(£^) (1) 

f 

where 6 Is the universal gravitational constant, given numerically by 
G S 6.6732 X 10"^^Nm\g"^. 

The force F experienced by P in the presence of P is 


F = - F (2) 

in conformity with Newton *s third law. 

Example ; If m is the mass of the earth, m the mass of the moon, and £ 

the vector from the mass center of the earth to the mass center of the moon, 

and if the numerical values of m, m and |£| are given approximately by 
m = 5.97 X lO^^kg, m * 7.34 x lO^^kg and |£j ® 3.844 x 10®m, what is the 
magnitude of the force ^ exerted on the moon by the earth? 

As Eq. (1) applies to particles, it can be used for the purpose at hand 

only with the supposition that the earth and moon can be replaced by 

The superscript 2 on a vector indicates scalar multiplication of the vector 
with itself, l.e., squaring the magnitude of the vector. 

^E.A. Hechtly, "The International System of Units: Physical Constants and 

Conversion Factors," NASA SP-7012, Revised, 1969. 
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2.1 


particles situated at their loass centers and having masses m and a, 
respectively. In that case, 

|F^| " Gmmlf^j (£/ ) 

( 1 ) 

S (3.98602 X 10^Wkg"^)(7.34 x lO^^kg) (3.844 x 10®m)’^ 

S 1.98 X 10^0 N (3) 

The validity of this supposition is examined in Sec. 2.4. 
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2.2 Force exerted on a body by a particle 

The system of gravitational forces exerted by a particle P of mss m 
on the particles of a (not necessarily rigid) body B is equivalent to a single 
force F vdiose line of action passes through P, but not necessarily through 
the mass center 5 of B. If B consists of particles P^, «..9 P^ of masses 


m«. 


. , m^, and if jg^ are the position vectors of Pj^, P^^ 


relative to P (see Fig. 2. 2.1), then F is given by 


N 


F - -Gm E 

i-1 ^ 


( 1 ) 


whereas, if B is a continuous distribution of matter, and p is the mass 
density of B at a generic point P of B, £ is the position vector of P 
relative to P (see Fig. 2.2.2), and dt is the length, area, or volume of a 
differential element of the figure (curve, surface, or solid) occupied by B, 
then F can be expressed as 


F - -Gmy^(£^)“^^^pdT 


< 2 ) 


Once the line of action of the force F in Eq. (1) or Eq. (2) has been 

1 

established, it is always possible to locate a point B on this line such 
that the force exerted by P on a particle whose mass is equal to that of B 

I 

and which is placed at B is equal to the force £ exerted on B by P. The 

f 

point B , called the center of gravity of B for the attracting particle P, 

t 

does not in general coincide with the mass center of B; and if R is the 

_ • 
distance from P to the center of gravity B ,, then 
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Figure 2.2.2 
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t 

R 


/ Gam 

\iii; 


(3) 


where P is obtained from Eq. (1) or Eq. (2)» and m is the mass of B. 

Derivations ; The force d£ exerted by F on a differential element of 
B at point P acts along the line joining P to P and is given by 

dP “ -Gmp(£^)”^^^pdT 

( 2 . 1 . 1 ) 

By definition, two systems of forces are equivalent if they have equal 
resultants and equal moments about one point. Now, on the one hand, the 
resultant of the system of forces exerted by P on all the differential 
elements of B is given by 


m -6m^£(£^)”^^^pdT 


and the moment of the system of forces about point P is equal to zero, because 
dP, acting along the line joining P to P, has zero moment about P. On the 
other hand, the resultant of a system of forces containing but one force F 
is £, Itself, and the moment of this system about point P is zero if the 
line of action of F passes through P. Hence, if F is given by Eq. (2) and 
acta along a line of action passing through P, then F is equivalent to the 
system of forces exerted by P on all differential elements of B. 

A parallel proof may be constructed for Eq. (1), which replaces Eq. (2) 
when B consists of a finite number of particles. 

The truth of the assertion that the line of action of the force F does 
not necessarily pass through B is most easily demonstrated by an example; 
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and Eq. (3) follows directly from the definition of R and Eq. (2.1.1). 

Example ; A uniform thin rod B of length L and mass m is subjected to 
the gravitational attraction of a particle P of mass m as shown in Fig. 
2.2.3. The system of forces exerted by P on the particles comprising B 
is to be replaced with a force, £ whose line of action passes through P> 

f _ 

and an expression is to be found for the distance R between P and the 

? 

center of gravity B of B for P. 

If the rod is regarded as matter distributed along a straight line 
segment 9 then the position vector £ of a generic point P of the rod relative 
to P can be expressed as 


£ = yb - ab 2 


where and b 2 are unit vectors directed as shown in Fig. 2.2.3 and where 
y varies from 0 to L. The mass density p of B at P is then equal to m/L, 
and a differential eletnent of B has a length dy. Hence the force F, 
resolved Into componants parallel to^^^ and b 2 , can be written 


— “ ^1—1 ^2—2 “ J ^^—1 ~ dy (5) 


and integration yields 


Fj^ - -Gmm ^1 + ^ | ^a(a^ + 


F 2 * Gmm/a(a^ + 


‘8 
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From Fig. 2.2.3 it can be seen that the distance e between the mass 

center B of B and the interscetion of B with the line of action of F is 
related to F- and F« by 



or 


e 

L 





(9) 


The ratio e/L from Eq. (9) is plotted vs. a/L in Fig. 2,2.4. Since two 
bodies cannot occupy the same point in space, the limiting case a/L ” 0 must 
be excluded from consideration. However, as this limit is approached, the 
ratio e/L approaches the value 1/2 and the line of action of F approaches 
coincidence with B. For any finite value of a/L, Fig. 2.2.4 indicates that 

— sfe 

the line of action £ cannot pass through both P and B . 

» _ I _ 

The distance R from P to the center of gravity B of B for P is given by 


R 


(3) 


Gmm 


1/2 




It is clear from this example that the location of the center of gravity 

I 

B is not in general a property of body B alone, but depends also on the 
position of the attracting particle P. 
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Figure 2.2.4 
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2.3 Force exerted on a small body by a particle 

When a body B Is subjected to the gravitational attraction of a particle 
P removed so far from the mass center B of B that the largest distance from 

ft 

B to any point P of B is considerably smaller than the distance R between 

mmmm 

P and B , a useful form of the expression for the force F given in Eq. (2*2 •2) 

can be found as follows: Replace £ with the sum of R, the position vector 

* — 

of B relative to P, and £, the position vector of P relative to B (see Fig. 
2.3.1), and then expand the integrand in ascending powers of lr|/R to obtain 



where f_'‘ ^ is a collection of terms of i degree in lr^|/R» ni and m are the 
masses of B and P, G is the universal gravitational constant, and is a 

— ft 

unit vector directed from P toward B , so that 


R = Ra^ (2) 

(2) 

In particular, ^ is given by 

— => {y [tr (^) - 5£j^*_I*aj^]aj^ + 3j[*aj^} (3) 


ft 

vdiere is the inertia dyadic of B for B , and tr(I) denotes a scalar in- 
variant of called the trace of 1 and defined in terms of any mutually 
orthogonal unit vectors n^^, 1 I 2 , and n^ as 

tr(l) A Si* I.’ El + E 2 *i’E 2 ^ — 3*i*— 3 
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Eq. (1) suggests two approximations which we shall find useful in the 
sequel: 


F 



Gmm 


( 5 ) 


and 


F = 




( 6 ) 


( 2 ) 

Two expressions for £ which sometimes furnish convenient alternatives 

to the one given in Eq. (3) are obtained as follows: Introduce unit vectors 

&2 and such that a^, & 2 , and a^ form a dextral, orthogonal set, and let 

— 1* —2* —3 unit vectors respectively parallel to principal axes of B 

* 

for B and also forming a dextral orthogonal set. Next, define and Ij 
as 



—j — ^ 


(j,k 1,2,3) 


and 


(j = 1,2,3) 

Finally, let 

^ij =1,2,3) 

( 2 ) 

Then ^ may be written either as 


( 7 ) 


( 8 ) 


0 ) 
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2,3 


:C2) _ 3 n 


_2 2 ^^22 ^33 ■ ^ 21-2 ^ 31 % 

laR ^ J 


( 10 ) 


or as 


:C2) „ _3_ 

' »2 
mR 


|I 1 j^( 1 - 3 gJ^) + I2CI-3CJ2) + ^-1 


+ Il3_C2j_Cj_j_ + 12 ^ 22^12 ^^ 23 ^ 13^-2 


■*■ ^^ 1 *^ 31^11 ^ 2 ^ 32^12 ^ 3 ^ 33 ^ 13^-3 


( 11 ) 


The relative simplicity of Eq. (10) is a result of the use of time-dependent 
moments and products of inertia tsee Eq. (7)]. By way of contrast, the 
principal moments of Inertia appearing in Eq. (11) are constants, and the 
orientation of B relative to a^, £ 2 * £3 comes into evidence through 

the direction cosines relating the two sets of unit vectors a^, ^ 2 * £3 

^1* -2' -3* 

Derivations ; Replacing £ in Eq. (2.2.2) by R + £ (see Fig. 2.3.1) 
provides 




F - -Gm / (R + r) • (R^ + 2R«r + r^)"^'^^pdT 


( 12 ) 


In terms of the normalized position vector £ defined by 


£ A r/R 


( 13 ) 


and the unit vector a^ satisfying 
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® R/R (14) 

( 2 ) 

the force becomes 

I “ - % /(ai + q) (1 + 2aTa + (15) 

(12,13,14) vrJ ^ ^ 

Application of the binomial series expansion 

(I4x)^ = 1 + nx + n(n-l)x^ + n(n-l) (n-2)x^ + . . . (16) 


(valid for x < 1) to the exponentiated quantity then yields, for 
2 

2 aj» 3 _+^ < 1, the expression 


F 

(15,16) 



[1- |(2a^*a. + /) 


+ |^(aj^*£)^ + ...] 


+ <[[1 - 3a^*3^ + ...]}pdT 


(17) 


where dots represent terms of degree three and higher in . This expression 
can be simplified by taking advantage of the fact that B is the mass center 
of B, for this means that 



0 


(18) 


so that 


F 

(17,18) 


(ka 


fi‘i - 1 


15 2 

yaiCaj^’a) - 3^[a^*a^]pdT + ... (19) 
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2.3 


Furthermore, replacing c[ with r/R and observing that 


f 


pdT = m 


one arrives at 


( 20 ) 


F 

(19,20) 


Gmia 


3Gm 

2R^ 


^aj^r^pdr - Sa^^aj^^Mi 


rrpdT'a, 



« • • 


( 21 ) 


The Integrals appearing in Eq. (2) can be related to inertia properties 

* 

of B by introducing two quantities: the inertia dyadic ^ of B for B , de- 

fined as 


I 4y*(Ur^ - 


rr)pdT 


( 22 ) 


where U denotes the unit dyadic, that Is, a dyadic such that, for any vector 
V, v*lJ = U*v = v; and the trace of I^, defined in Eq. (4) . For it follows 
from Eqs. (22) and (4) that 


tr(l) 



+ n„»rr»n^ 




2 #r pdT 


nj*rr*nj)]pdT 


(23) 


so that 


f 


^ A 1 

r^ pdT = — 

(23) ^ 


trCD 


(24) 
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* fl tr(l) 

rrpdT = U I r pdT - ^ = U — _I 

C22) J (24) ^ 


Consequently 


t, - - ^ £l - ^ Itt® - 5a,-I-.Ja, - ^ I-a, + ... (2«) 

(21,24,25) r ^ 2R^ ^4 _ -i 


and the equivalence of Eqs* (1) and (26) becomes apparent after using £q« (3) 
and recognizing that the three dots in Eq, (26) stand for terms of degree 
three and higher in jr^j/R, these terms being represented in Eq. (1) by 
for i ^ 3 . 

Referring to Eq. (7), one can express as 


- " ^12-1-2 * ^13%-3 


^ 21 — 2 % ^ 22 — 2—2 ^ ^ 23 ^ 2 % 


+ '^ 21 — 3^1 ^ 32%^2 ^ 33%^3 


from which it follows that 


" ^11% ^21-2 ^31-3 


Substituting from this equation into Eq. (3) and noting that 


trCD 


(4,7) 


I + I + I 
11 22 ^ 33 
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^ile 


-1--1 ^11 
one arrives at Eq. (10). 

Finally, Eq. (11) may be obtained from Eq. (10) after observing that, 
in view of Eq. ( 8 ) and of the assumption that b^, ^2* ^3 parallel 

to principal axes of inertia of B for B , I can be expressed as 

I = + hhh * ^^3 


so that 


•■ll ^jy ^l%‘-l-l’-l 


+ I^a- •'b,^bA*a- + l^a-*b^b«*a- 
2—1 — 2—2 —1 3^1 «=- 3=-3 — 1 


(9) 


Vll h^l2 ^3^13 


and, similarly. 


^21 “ ^ 1 ^ 21*^11 * ^ 2 ^ 22^12 ^ 3 ^ 23^13 

and so forth. Substitution into Eq. (10) then yields Eq. (11) if one makes 

2 2 2 

use of the fact that + C ^2 ^ i ■ 1» 2, 3. 

Example ; An approximate expression is required for the force F exerted 
on a uniform thin rod B of mass m and length L by a particle F of mass m 
located relative to B as shown in Fig. 2.3.2, with R » L. Body B is to be 


19 
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idealized as matter distributed uniformly along a straight line 
segment. 

When i etc, are omitted, Eq. (6) provides the required ap- 

proximate force expression. If Ib^, and are Introduced as a set of 

dextral, orthogonal unit vectors, with parallel to the rod axis, then 

A 

the inertia dyadic ]t of B for B is given by 

i ■ + W 

which, when substituted into Eq. (3), introduces into the force expression 
the dot products — 3*— 1* evaliuite these, let ^ be the angle 

between and a^^, and note that only these two unit vectors have physical 
significance inherent to the problem, so that b 2 > b^, a 2 and a^ are defined 
arbitrarily, as shown in Fig, 2,3. 2, It follows that 


1 * ^ 

(28) 

^3*^1 “ ° 

(29) 


and 


sin^ ^ 




so that, after using 


^2 “ ®2 ' I ' 


(30) 


(31) 


21 
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and omitting \ etc., one obtains 


^ _ WWl 

SB 5* W mm 

“ (6) (6,30,31) R' 


— { r 2 I 2 

^ a. 1 + ^ (2-3sin^f) - a. —• sin2»{) (32) 

R I 8R ^ 8R 


This result could have been obtained as easily from Eqs. (1), (.6)', and 
(11), with the substitution of directlW cosines available from Fig. 2.3.2 as 


cosf -slni|) 0 
slntj^ co8tj> 0 
0 0 1 


11 

"l2 

"l3‘ 

21 

•^22 

^23 

31 

^32 

^33 . 

I. I. 


_ oL^ 

'•I'rt 


* — 

When the rod Is aligned with the line joining Its mass center B to P, 
so that a 0, the force becomes 


Gmm In, L \ 
(32) R^V 4R^)“^ 


and when the rod Is perpendicular to this line, so that t{) = ir/2 rad., 

/ .,2 V 


(32) R 




Since in these special cases lines passing through B and paralleling a^, 
^2> ^1^8 a^ ere principal axes of B for B , Eqs. (33) and (34) could have 
been obtained most directly by using Eq. (10), with = l^j^ °° 0 and for 

t - 0, 

^11 “ ^22 “ ^33 “ ^ 


22 



2.3 


whereas, for i|» >■ tr/2. 


I 


11 



^22 " ^33 



Note that the term (GmmL^/8R^)sln2^£2 (32) represents a force 

directed normal to the line joining P to B . Such force components cause 

bodies to move in orbits differing from the Keplerlan orbits associated 

with particles ; but this effect is so small that it can generally be Ignored 

even in high precision orbital calculations. 

Certain differences between the Example in Sec. 2 .2 and the present 

example should be noted carefully. For instance, Eq. (32) is valid only 
L 

when ^ "~2 3 , as required for the underlying series expansion 

4R 

in Eq. (17), whereas Eqs . (2.2.5)-(2.2.7) in the Example in Sec. 2.2 are 
not subject to any such restriction; and the line through P and normal to 
B passes through the end of B in Fig. 2.2.3, but need not do so in Fig . 

2.3.2. For the special case in which P in Fig. 2.3.2 lies on the normal 
to B passing through the end of B, as in Fig. 2.2.3, cos ^ and sin 4* are 
given by 

cos tp * L/2R (35) 

and 

sin 4» - [1-(L/2R^)3^^^ (36) 


but the form of Eq. (32) applicable to this case would not be obtained by 
simply substituting these expressions into Eq. (32); rather, it is necessary 


23 


2.3 


to expand sin ij; in ascending powers of L/R and then to drop terns of degree 
three or higher after substituting into Eq. (32), which gives 


F 



(37) 
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2.4 Force exerted on a small body by a small body 

* — * 

When the distance R between the mass centers B and B of two (not 
necessarily rigid) bodies B and B substantially eiKceeds the greatest dis- 
tance in either body from the mass center to any point of the body, the 
system of gravitational forces exerted on B by B has a resultant £ which 
can be expressed as 


F 



00 00 . 00 

+ 1: E 

i-2 i-2 J-2 / 


( 1 ) 


where a^ is a unit vector directed from B toward B , 6 is the universal 
gravitational constant, and m and m are the masses of B and B, respectively, 
and where is a collection of terms of i^^ degree in |r^|/R, is a 

collection of terms of 1^^ degree in lic|/R» and is a collection of 

terms in the product (|xl/R)^( |r ] /R)^ , with _r and r respectively position 
vectors of generic points of B and B relative to B and B . In particular,. 



{tr(^)-5aj^*j[‘aj^]aj^ + 31 


( 2 ) 


and 



3 

2 


I tr (I) -Saj^ *l*a,j^] 


+ 


3I*aj^ 


(3) 


* — 

where and ^ are the inertia dyadics of B for B and of B for B , 

respectively. 

If a ^2 defined so as to establish a dextral, orthogonal set of 

( 2 ) —( 2 ) 

unit vectors a^, £ 2 , a^, then f^ and can be expressed in terms of 

these unit vectors and the moments and products of Inertia of B and B 
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for axes paralleling a^, a^, and and passing through the laass centers 

of the individual bodies. To this end, I,, and I., are defined as 

3 K 3 * 



4 


(j,k - 1,2,3) 


( 4 ) 


and 


a.k- 1,2,3) <5) 

^23 ^23 

after which JE'^ and ' nay be written 

^ ft <*22 + *33 - **ll)il + *21*2 * *3li3] <« 

and 

^ [2 ^^22 ■*“ ^33 ■ ^^11^% ^21-2 ^31-3] 

( 2 ) ( 2 ) 

Alternatively, f^' and * can be expressed in terms of principal 

^ >•«» - 'it 

moments of inertia of B for B and of B fox B , To accomplish this, two 

sets of dextral, orthogonal unit vectors, bj^, b 2 » b^ and bj^, b^, para- 

^ — 

lleling principal axes of B for B and of B for B , respectively, are 


introduced, and 1^, I^, 

'13 • '13 

defined as 


*3 ■ 

(3 - 

1,2,3) 

(8) 

*3 

(j “ 

1,2,3) 

(9) 
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(i,4 ■> 1,2,3) 


'ij ^ 


(i,j - 1,2,3) 


Thus one obtains 


— ” 2 2 + l2<1^3Cj^2) ^3 ^^"^^13^^% 

mR 


^ ^^ i *^ 21^11 ^ ^ 2 ^ 22^12 ■*“ ^ 3 ^ 23 ^ 13^-2 


■*■ ^^ 1 ^ 31^11 ■•■ ^ 2 ^ 32^12 ^ 3 ^ 33 ^ 13 ^% 


“ z\ i| IIiCl-3cJj_) + I 2 (1-30^2^ + 

mR 


'*' ^^1^11^22 ^ ^2^22^12 '*' ^3^23^13^-2 


+ IIxC3j^C^2 + ^2*^32^12 ■*“ ^3^33^13^-3 


A useful approximation of j? In Eq. (1) may be obtained by defining £ 


such that 


P S i A - (a- + f + f ^^^) 
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with and substituted from either Eqs. (2) and (3) , or Eqs. (6) 

and (7), or Eqs. (12) and (13). 

Derivations ; To establish the validity of Eqs. (1), (2), (3)* (6) , (7), 
(12), and (13), expressions given in Sec. 2.3 may be used to represent the 
force exerted on B by a differential element of B located at a generic 
point P of B (see Fig. 2.4.1), and the total force £ applied to B by B 
can then be obtained by Integrating over the figure of B. Specifically, 
if R is the distance between F and B , is a unit vector directed from 
P toward B*, p is the mass density of B at P, and dT is the volume of a 
differential element of B at P, then, from Eqs. (2.3.1)-‘(2.3.3) with R, 
and m replaced by R, a^, and pdt, respectively. 


Gm/- 

:r2\% 


I Itr(l) - Saj^* 


■Si'ii 


+ 3I*8j^ 


where the three dots represent terms of third or higher degree in |r_|/R. 

A ~ 

Now, if ^ denotes the position vector of B relative to P, and R the 

^ A 

position vector of B relative to B , then (see Fig. 2.4.1) 





Hence 



( 16 ) 


F » Gm /(R + r) (R^ + 2R*r + r^)"^^^ pdx 
(15,16) J 


- Itr(I) - 5aj^‘X*a^]aj^ + + ..^pdT (17) 
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The first Integral in this equation has precisely the same form as the 
integral in Eq. (2.3.12). Consequently, proceeding as in that derivation, 
one obtains results analogous to Eqs. (2.3.26), namely 


6m 


nj ' (R + r)(R^ + 2R*£ + pdt 


Gmm / . 1 




i> 


(18) 


where the three dots represent terms of third or higher degree in Ijtj/R. In 
the second Integral in Eq. (17), R and may be replaced with R and 
respectively, because every term in the ^integrand involves quantities of 
second or higher degree in |r_| , so that no terms of interest for the purposes 
at hand are lost through this replacement; and, once the replacement has been 
made, the portion of the integral displayed explicitly in Eq. (17) can be 
evaluated readily. Thus, one obtains 


F - 
(17,18) 




j Itr(l) - 5aj^*Z*aj^ia.j^ 


+ 3l*aj^ + 


• • • 


mR 


> 


(19) 


where the three dots now represent terms of third or higher degree either in 
jr^l/R or in |r|/R or terms Involving the product (|£|/R)^ with 

neither i nor j equal to unity, because £ and £ are drawn from the mass 
centers of B and B. Eq. (1) now follows directly from Eq. (19) together 
with the definitions in Eqs. (2) and (3). 
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The relationship between Eqs.(2), (6) , and (12) is completely analogous 
to that between Eqs. (2.3.3), •—(2.3.10), and (2.3.11); similarly, for Eqs. 

(3), (7), and (13). Hence, to establish the validity of Eqs. (6), (7), and 
(13), one can proceed exactly as in the corresponding derivations in 
Sec. 2.3. 

Example ; An approximate expression is required for the force exerted 
by a homogeneous, oblate spheroid B on a homogeneous rectangular paralleplped 
B having the dimensions shown in Fig. 2.4.2. The approximation of £ denoted 
by F in Eq. (14) is to be employed, and the three additive terms in this 
expression are to be compared with each other. For purposes of numerical 
comparison, the values 

g 

a “ 6*38 X 10 m 

e = 6.36 X 10® m 

a = 16 m 
3 < a 

R « 6.6a - 4.2108 x lO^ m 

are to be used. The system then roughly approximates that of the earth (B) 
and a large, synchronous-altitude artificial satellite (B). 

Eqs. (12) and (13) provide a convenient point of departure for tl^ 
required comparisons. The principal moments of inertia appearing in these 
equations are given by 


A satellite in circular orbit at synchronous altitude orbits the earth 
every twenty-four hours, so that when in ah equatorial orbit the satellite 
remains above a fixed point on earth. 
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m ^ + ot^) (20) 

^ 2 - 13-12 

for the rectangular parallelepiped B, and by 

Tl - I (a^ + a^) (22) 

I 2 “ I 3 “ f (a^ + 6 ^) <23) 


for the spheroid B. Later work is facilitated by eliminating B and B through 
the introduction of the eccentricity e of the spheroid B with the substitution 

B^ = (1-e^) (24) 


and a similar quantity e for body B with the substitution 


B^ = (1-e^) 


(25) 


( 2 ) 

The vector f then becomes 


-'\ih ^ ^ 


+ 2[e23Cj^j^ + f^22^12 ^23^13 “ 2 ^^22^12 ^23*^13^^-2 


+ 2{C23_Cj_j_ + C32Cj_2 + ^22pl3 ~ 2 ^^32^12 ^33^13^ ^-3^ 


(26) 
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org after simplification by means of Eq. Cl»2«14) j 


:C2) 


(26) 4E 


2 2 

. 2 [2 * ^21^11-2 ^31^11%] 


( 27 ) 


By using Eq. (1.2.15), one can now express the magnitude jf' '| of f_'‘ ^ as 


|£<«l . a!|i (1 - 20?, + 

(27) 8R^ 


(28) 


Similarly, 



3a^e^ 

5R^ 


(1 


^^ 11>^1 


^ 21 *^ 11-2 


Sl^ll% 


(29) 


and 

lj(2)| ^ .. S (1 _ 2C^ + 5^,)^^^ (30) 

lOR*^ 

2 4 1/2 

Fig. 2 .4 .3, showing a plot of the function (1 - 2x 4* 5x ) versus 
X in the range -1 ^ x < 1 , can be used to find extremals of |f ^ | and 

I substituting and for x. By reference to Eqs. (28) and (30) 
it can thus be concluded that 


2^2 

a e 


r(2)| 


4v3r' 


2 - 


2,-2 

(X e 

4R“^ 


and 


ii'"’ 


— 2—2 
. 3a^g 

— 2 
5R"^ 


(31) 


(32) 
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and substitution of the given numerical values then yields 


1,(2) 

|J(2) 


max 

min 


Ml 

12 - 2^2 


; lO'V 

12(6,38)^xl0^2<.j006> 


< 10 


-9 


(33) 


and 



max 



iSS 


5R‘ 




3 (.006) 
5(6.6)^ 



(34) 


Eq. (34) shows that, for a large, synchronous-taltltude satellite of the 

earth, the leading term in the series for F. in Eq. (1) has a magnitude far 
—( 2 ) 

exceeding that of ^ > a vector which reflects the oblateness of the earth 

through the presence of e in Eq. (30) ; and (see Eq. (33)) | f.^^^ | is much 
smaller than even j * \ . 

( 2 ) 

It is noteworthy that the tiny term ^ ' (see Eq. (27)) depends on 

^11’ ^21’ ^31’ hence on the orientation of the satellite relative to 

—( 2 ) — — — 

— 1’ — 1 * —3* X (®®® (29)) involving €2^^, and 

depends on the orientation of the satellite's orbital plane relative to 

the spheroid. For an equatorial orbit, reduces to 


J(2) 


3a2r2 


(29) 


5R' 


2 —1 


and the addition of this vector to the term in Eq. (14) is of little 
consequence (see Eq. (34)). By way of contrast, if B moves in a non- 
equatorl^ orbit, oblateness effects may be significant despite the fact 
that is small in comparison with because then has components 

perpendicular to 
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2.5 Centrobaric bodies 

I 

As noted in Sec. 2.2, the center of gravity B of a body B for an 
attracting particle F does not In general coincide with the loass center 
B of B. However, there exist bodies for 'tdiich the center of gravity and 
center of mass necessarily, coincide. Such bodies are called centrobaric. 
Thus a body B of mass m Is centrobaric if the force j? exerted on B by 
every particle P of mass m is given by 


F 


Gmm 


.2 -1 


( 1 ) 


where 6 is the universal gravitational constant, R is the distance between 

— » ik “ft 

P and B , and a^ is a unit; vector directed from F toward B . 

Centrobaric bodies cah be fodnd in a variety of interesting shapes and 

mass distributions (see Prob. 2f ) ; but all such bodies possess the followlne 

property: A centrobaric body has the seme moment of inertia about every 

line passing through Its mass center. In other words, the inertia ellipsoid 

of a centrobaric body for the mass center of the body is a sphere. However, 

not every body possessing this property Is centrobaric (see Prob . 2e) . 

Derivations : The existence of centrobaric bodies is most easily proven 

by citing a specific case (see the Example, in which a solid sphere is 

shown to be centrobaric if its mass density at any point is a function only 

of the distance from the point to the center of the sphere) . 

To prove that a centrobaric body B has the same moment of inertia about 

* 

every line passing through its mass center B , it is sufficient to show 

In the classical literature, the term center of gravity is sometimes so de‘<^ 
fined that only a centrobaric body has a center of gravity; then the center 
of mass coincides with the center of gravity whenever the latter exists. 
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that the quantity has a value independent of where X is the 

* 

inertia dyadic of B for B « Nov, taken in conjunction, Eq* (1) and the 
series expansion for F given by Eq. (2«3.1) imply 



i=2 


The terms in this summation are independent of each other in the sense 

(2) ~2 ( 3 ) -3 

that _f is proportional to R , f is proportional to K , and so forth, 

and R may be taken to be arbitrarily large (see Eq. (2.3.3). It follows 

that they vanish separately, that is, that 

=0, i * 2 op 


Hence 




( 2 ) 

and, using the expression for f^ given in Eq. (2.3.3), one finds that 


- 3 


Since tr(^) is an invariant, the moment of inertia ® about the 

A 

line that passes through B and is parallel to has been shown to have the 
same value for all orientations of a^ relative to B. 

Example ; To show that a solid sphere S of mass m is centrobarie if 
the mass density p at a point P depends only on the distance r between P 
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and the center S of S, the force F exerted on S by a particle P of mass m 
will be calculated. 

The position vector p of P relative to P, expressed In terms of the 

«— ic 

spherical polar coordinates r, 6, ip, the distance R between P and S , and 
the unit vectors a^, a^, shown In Fig. 2.5.1, is given by 

=* (R + r cos ip)a^ + X sin ip sin 6 £2 ^ ^ ® fi.3 

The volume dr of a differential element of S is 

dT " r^ sin ip dr d0 dip (3) 


and with p = p(r) the mass m Of S can be expressed as 


r 2 f 

m “ I pdT “ / X p I sin ip I d0 dip dr 

J •'0 


2ir 


or 


m * Air / 


r^pdr 


where a Is the radius of S. 

The force F exerted by P on S is given by 


ifi 


- Gm I jg*(£^)“^^^pdT 


C2.2.2) 

or, In somewhat more explicit terms, by 


(A) 


C5) 
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Figure 2.5.1 
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F 

(2,3, 


- Gm r^p f sin4(R^+2Rr cos4+r^)"^^^ /* £ d6 dif) dr 

5) Jo «^0 


( 6 ) 


Using Bq. (2), one obtains for this inner integral in Eq. (6) 


/ 2‘ir 


I(R+r cosi{»)aj^ + r8in^8in0a2 + rsini^fcosSa^WQ ~ 2iT(R+rcosij^)aj^ (7) 


so that 


F “ -2'n’Gm f r^p f sin^(RV2Rrcosij>+r^)”^^^(R+rcost|))dt[)ir a. 

(6,7) -'0 J) ^ 


( 8 ) 


Integration is now facilitated by the introduction of a new variable v 
defined such that 


A R^ + 2Rr coB\p + r^ 


(9) 


When this definition and the implied relationships 


sin^di)> ” - vdv/(Rr) 

(9) 


( 10 ) 


and 


R + r cost|» - (v^ + R^ - r^)/(2R) 
(9) 


( 11 ) 


are substituted into Iq. (8), there results 
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_ /-a 

irGm 

•'0 


/ ® /•(R-r), 2^„2 2v 

rp f dv dr a. 

•'(R+r) 


AttGui / ® 2 , Gm m 

— 5“ / r pdr a, - ” —5- a^. 

R^ •'0 ^ (4) R^ ^ 


( 12 ) 
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2.6 Moment exerted on a body by a particle 

The system of gravitational forces exerted on a (not necessarily rigid) 
body B of mass m by a particle P of mass m produces a moment M about the 
mass center B of B, given by 

M = - R X F (1) 


where F Is the resultant of the system of forces [see Eqs. (2.2.1), (2.2.2), 

^ 

and (2.3.1)J and R is the position vector of B relative to P. If the 

— k 

distance R between P and B substantially exceeds the greatest distance 

k 

from B to any point F of B, this moment can be expressed as 


M 


3Gm 


3 + 


Gmm 

R 


O® > V 

1-3 


( 2 ) 


where Is a unit vector directed from P toward B , G is the universal 

k 

gravitational constant, Is the inertia dyadic of B for B , and the dimen- 
sionless vector m^^^ is a collection of terms of 1^^ degree In |r.| /R, with 

k 

the position vector of a typical point P of B relative to B . 

Eq. (2) suggests the approximation 


M ~ M A 

- -L 1 


(3) 


Expressed in terms of scalars ” 1>2,3, defined as 

(j,k - 1,2,3) (4) 

the vector M is given by 
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M 


3Gm 


<^ 21 % 


^ 31 - 2 ^ 


(5) 


Alternatively^ one can introduce a deKtral set of orthogonal unit vectors 

62 , ^2 parallel to the principal axes of inertia of B for B , and 

express H in terms of the principal moments of inertia I 2 , of B 
* 


for B and the direction cosines C 


ij’ 




1,2,3 defined respectively as 


(j - 1,2,3) ( 6 ) 

and 

I ’ 

(i,j - 1,2,3) (7) 


thus obtaining 


3Gm 


M * 3 t^l^V^ 2 ^^ 12^13 - 2 ^^r^ 3 ^^ 13 ®ll - 3 ^^ 2 "^ 1 ^^ 11 ^ 12 ^ 


( 8 ) 


(5), despite its apparent simplicity, is less useful than Eq. ( 8 ), 
both because the products of inertia l 2 j^ and vary with the orientation 
of relative to B and because the rotational equations of dynamics are 
generally most easily fomulated in terms of the vector basis jb 2 » b^. 

It should be noted that M, vAiich by Eq. (3) is the first term in Eq. 
(2), can vanish when M does not vanish. Hence, when using M as an approxi- 
mation to M one is not necessarily retaining the largest term in the series 
expansion for H. Specifically, M vanishes whenever parallels a principal 
axis of inertia of B for B , but this state of affairs need not produce a 
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zero value for M (see Prob. 2g) . Moreover, H is identically zero for any 
body with spherical Inertia ellipsoid, whereas M is identically zero only 
if that body is also centrobarlc (see Sec. 2.5 cmd Prob. 2e) . 

Derivations ; If B is a continuous distribution of matter, then the 
resultant F of the gravitational forces exerted on B by P can be expressed 
in terms of the force dF exerted by F on a differential element of B at a 
generic point P of B as 


F 



(9) 


and, if r is the position vector of P relative to the mass center B of B, 

A 

then the moment M of the system of gravitational forces about B is 


M 




( 10 ) 


Substituting 


r - - R + £ 


( 11 ) 


where R is the position vector of B relative to P and £ is the position 
vector of P relative to P (see Fig. 2.6.1), yields 


M - / (- R+£)x^ - -RXF+/£X^ 

(10,11) J (9) J 


( 12 ) 


But ^ parallels £, so that £ x ^ > 0 and 


M - - R X F 

( 12 ) 


( 13 ) 
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confirming Eq. CD* Alternatively, this result can also be inferred fromi 
Sec. 2.2 and from the fact that two equivalent systems of forces have equal 
moments about every point. Hence It applies also ^en B consists of a | 
finite number of particles, Pj^, ..., I 

Eq. (2) for H is to be proven under the s^ne restriction^ on R ai^ c|n 
the dimensions of B as apply to the expression for F in Eq. (2.3.1). Uncfer 

i ] '/) 

these circumstances, substitution from Eqs. (2.3.1) and (2.3.2) into Eq. | 

(13) yields 

M - ^ a, X f; £<« - (a, X + f; a, X f<«) (14) 

^ 1-2 ~ R -1 - -1 - 

and use of Eqs. (2.3.3) and (2.3.4) then gives 


„ 3Gm „ , . Gm m 

M “ — 5 — a, X I»a, + 


r3 H i.3 


T. a, X f 


( 1 ) 


(15) 


which is equivalent to Eq. (2) if m^^^ is defined as 




9 • • • > 


(16) 


Verification of the equivalence of Eqs. (5) and (8) with Eq. (3) require^ 
the substitution 


—1 ” ^ii^ 

^ ( 1 . 2 . 2 ) 



After setting 


I 


(4) 


ij-l=j 
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one then atrlvee at Eq. (5) > and using 


I 


( 6 ) 


+ ^ 2 ^ 21^2 + ^ 3 ^ 3 ^ 


( 19 ) 


one obtains Eq. ( 8 ) . 

Example ; The vector M defined in Eq. (3) is to be used to approximate 

— — * 

the moment H exerted by a particle P of mass m about the mass center B of 
a homogeneous, right circular cylinder B having moments of inertia J and I, 
respectively, about the sysmietry axis and about any line through B and 
normal to the S 3 nnmetry axis. The result is to be expressed in terms of two 
angles, and 6 , used to specify the orientation of the symmetry axis of B 
relative to a dextral, qrthogonal set of unit vectors S. 3 * Specifi- 

cally, as shown in Fig. 2.6.2, 6 is the angle between a^ and the symmetry 
axis of B, and (|) is the angle between intersection of the plane P 

passing through B and normal to a^ with the plane Q determined by the 

'ii 

symmetry axis of B and a line passing throu^ B and parallel to a^. 

Of the alternative expressions for M given in Eqs. (3), (5) and ( 8 ), 
the last is the most convenient, for suitable unit vectors b^, ^ 2 * ^3 

can be introduced readily, for example as shown in Fig. 2.6.2, where 
is normal to plane Q. Using Eq. (7), one then finds that 

= cos <j) 

C^2 “ “ ® si® ^ 

“ sin 9 sin <j> 
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and, in accordance with Iq. (6), 


Hence 


M “ (I - J) sin 8 sin 4>[b- cos 6 sin ()) + b„ cos ((>3 

“ ( 8 ) ^ 

If I ^ J, M vanishes only when at least one of the following conditions 
is fulfilled: sin <j> “ 0, sin 0 « 0, or cos 0 ® cos 4 “ 0. In all of these 
cases, the symmetry axis of B is either normal or parallel to the litte joining 

-—'k 

P and B , and considerations of symmetry indicate that M also vanishes under 
these circumstances. Thus it appears that M » 0 can imply M 0 for a 
particular body, although this implication is not valid in general. 
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2.7 Moment exerted on a small body by a small body 

When the distance R between the mss centers B and B of two (not 
necessarily rigid) bodies B and B substantially exceeds the greatest dis-- 
tance from the mass center of either body to any point of that body^ the 
system of gravitational forces exerted on B by B produces a moment M about 
B which can be expressed as 


^ X I.., + ^ E t 

^ i_3 » i-2 j-2 


( 1 ) 


where Is a unit vector directed from B toward B , ^ Is the Inertia dyadic 

— — 
of B for B , m and m are respectively the masses of B and B, 6 Is the uni- 
versal gravitational constant, and the dimensionless vector m^^^ Is a 

th 

collection of terms of l''” degree In |£|/R, while the dimensionless vector 
m^^^^ Is a collection of terms In the product (|rl/R)^(|r | /R)^ , with £ and 

£ respectively the position vectors of generic points of B and B relative to 

* -* 

B and B . 

The similarity of Eqs, (1) and (2.6.2) suggests that an approximate 
relationship similar to Eq. (2.6.3), namely 


R 


( 2 ) 


may prove useful. The vector M thus defined can be expressed as 


rj 3Gm _ >. 

- ” _3 ^^21% " ^31-2^ 

R 


( 3 ) 


or as 
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M - ^3“ ^12^13 -2^^r^3^^13^11 ■*■ -3^^2~^1^‘^11^12^ 

R 


(4) 


where I., I.., and G.. have the same meaning as in Sec. 2.6, and M, 

X 1 ij ij 

that Is, the first term in Eq. (1) , can be dominated by other terms In the 

series. A significant difference between Eq. (2.6.2) and Eq. (1) is that 

the former can be replaced with Eq. (2.6.1), whereas the latter cannot be so 

replaced, even if R and ^ are re-defined respectively as the position vector 

of B relative to B and as the resultant force exerted on B by B. 

Derivations : The validity of Eq. (1) can be established by using Eq. 

* 1 

(2.6.2) to represent the moment exerted on B about B by a differential 
eles^ent of B located at a generic point F of B (see Fig. 2*7.1) and then 

integrating over the figure of B. Specifically, if R is the distance 

^ ^ 

between P and B , a^ is a unit vector directed from P toward B , p is the 

mass density of B at P, and dr is the volume of a differential element of B 

at P, then, from Eq. (2.6.2), with R, and m replaced by R, a^, and 

pdT, respectively. 


H 



X I.*a^ + 


* • « 



(5) 


where the three dots represent terms of third or higher degree in |r^|/R. 
Because ^ consists of terms of second degree in |jr|/R, the integration in 
Eq. (5) can produce only terms of second and higher degree in |r^|/R, -and 
the substitution of a^ for and R for R cannot result in the loss of any 
terms of Interest. These substitutions permit the first term in parentheses 
to be removed from the Integrand, 'leaving as a factor an Integral equal to 
the mass m of B, confirming the first term in Eq. (1) ; and they yield the 
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first series in Eq. (1). The second series in Eq. (1) reflects the deviations 
of from and R from R. Every term in this series involves and 

it does not contain terms |.inear in |£|/R because r is drawn from the mass 
center of B. 

Eqs. (3) and (4) can be obtained from Eq. (2) by a procedure analogous 
to that used to derive Eqs. (2.5.5)' and (2.5.8) from Eq, (2.5.3). 

Example ! In the Example of Section 2.4 (see Fig. 2.4.2), an approxi- 
mate expression was developed for the force exerted by a homogeneous oblate 
spheroid B on a homogeneous rectangular parallelepiped B. Now an approxi- 
mate expression is to be obtained for M, the moment exerted by B on B about 
the mass center B of B. To this end, ap given by Eq. (4) is to be 
formed, and the magnitude of M is to be determined for 

m “ 6 X 10^^ kg. 

m •• 160 kg. 

a « 16 m. 

3 ■ 4 m. 

R ■ 4.2108 X 10^ m. 

G - 6.6732 X 10"^^ Nm^ kg“^ 

As in the Example of Sec. 2.4, the system roughly approximates the 
earth B emd a large, synchronous altitude artificial satellite B., 

The symmetry of B permits the substitution of I 2 for in Eq. (4), 
furnishing 


M m 
- 


« 1 - 


^ 2 ^‘^ 11 ^^ 13^2 ■ 


"= 12 ^ 3 > 


( 6 ) 
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from which 


lal ■ ^ <h-h>hAyi2>""‘ 


<7) 


Variations in the magnitude of M resulting from changes in the relative 
orientation of B and B thus depend only on changes in the angle between 
and and |k| attains its maximum value when this angle is equal to ir/4 rad., 
since the derivative of |h| with respect to vanishes when is 
Thus, for any body B with “ l 2 » 


lii 


max 


3Gm 

2R^ 


(^1 - 


V 


( 8 ) 


and for the given rectangular parallelepiped, for which 



( 9 ) 


and 


*2 ■ *3 ■ l2 

IhI becomes 
'—’max 

GSn(ot^-B^) ]im(a^~B^) 

3 “ 3 

or, after using tke given numerical values. 


|m| 

' — 'max 


ao) 


(11) 
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~ 6 .67 X 10~^^ X 6 X 10*^** X 160 X (256 - 16) j. 
8 X (4.2108)^ X 10^^ 


s24 


S 2.6 X 10 


-5 
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2.8 Proximate bodies 

When dealing with the gravitational Interaction of two bodies which are 
in close proximity, one can at times make effective use of results which, 
at first glance, may appear to be Inapplicable because 'they were encountered 
initially in connection with the analysis of the gravitational Interaction of 
widely separated bodies. For Instance, consider Eq. (2.7.2), which furnishes 
an approximation that becomes ever better as the distance R between B and 
B (see Fig. 2.8.1a) grows in comparison with the largest dimension of either 
body. Eq. (2.7.2) can be useful also when B and B are near each other (see 
Fig. 2.8.1b), provided B does not differ too much from a centrobarlc body 
(see Sec. 2.5) and R is sufficiently large In comparison with the largest 
dimension of B. For B then acts nearly like a particle of mass m situated 
at the mass center B of B; the resultant moment about B of the forces 
exerted on B by such a particle Is given to a good approximation by Eq. 
(2.6.3); and this equation is Identical with Eq. (2.7.2). 

Example ! Fig. 2.8.2 shows a particle P of mass m situated on the 
axis of revolution of a uniform oblate spheroid B of mass m, at a distance 
R from the mass center B of B. To explore the utility of Eq. (2.3.6) in 
situations Involving proximity of a particle and a body, the ratio F/F is 
to be plotted versus 3 for various values of e, where F Is the magnitude of 
£ as given in Eq. (2.3.6), F Is the magnitude of F^ as given In Eq. (2.2.2), 
and 3 and e are defined as (see Fig. 2.8.2 for a and b) 


3 4 b/R 

(1) 

£ 4 n - Cb/R)^J^/^ 

(2) 
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Figure 2.8.1 
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If and a^ unit vectors as shown In Fig. 2.8.2, and ** ^ 

then the associated moments and products of inertia of B for B are given by 


I I .1 

5 ’ ^22 ^3 


I (a^ + b^) 


^12 * ^23 " hi ■ 0 


Hence 


f(2) 

C2.3.10) 



/ 2 . 2 . 

(a - b ) aj^ 


and 


C2.3.6) 


( 1 , 2 ) 



F 4 lil 





- iiV 
5 l-e2 


(3) 


For the evaluation of the integral in Eq. (2.2.2), it is convenient 
to introduce the coordinates r, 6, and z shown in Fig. 2.8.3 and to note 
that £, the position vector of P relative to F, then can be expressed as 


£ ■ (R - z)aj^ + r sin 6 £2 ® S3 


while 


dT « r dd dr dz 
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As for p, the mass density of B, this is given by 


P 


3m 

47Tba^ 


Hence 


36mm 




(R-z)aj^+r8in0a2+rco86a2 


(2.2.2) 4iTba^ I(R-z)V3^^^ 


r dr dz d6 


where 


®1 " ®2 “ *1 " ^2 “ ^2 “ all-Cz/b)^]^^^ 


and» after carrying out the indicated integrations and using (1) and (2) to 
eliminate b/R and a/b, one arrives at 


F A |f| 


3Gmffi(l-e^) 


IH(1) - H(-l)1 

2$e 


(4) 


where H is defined as 


H(5) A sin"^ 


l-e^a-6€) 


Consequently, 


£ 

F 


(3,4) 3(l-e^) 


1 - 


3 


^ 1-e^ 


_ (1-g^) 


1/2 


2$e 


IH(1) - H(-l)3 


(5) 


62 



2.8 


In Fig. 2.8.4, F/F as given by Eq. (5) is plotted versus B for four 
values of e. Note that small values of 6 correspond to placing P at a great 
distance from B (see Eq. (1)). Hence it is not surprising that, regardless 
of the values of e, F/F approaches unity as 0 approaches zero. Conversely, 
values of 0 near unity represent situations in which the particle comes 
relatively close to the spheroid; and the error one then makes when using 
F in place of F can be seen to depend on e, larger departures of F/F from 
unity being associated with larger values of e, that is, with a more pro- 
nounced flattening of the ellipsoid (see Eq. (2)). However, it appears that 
even for 0 ~ 1.0, that is, for a particle that is nearly in contact with the 

spheroid, F/F departs only very slightly from unity so long as e < 0.10, 

* * 

which is the case, for example, when e •• e , where £ is the eccentricity 
of a spheroid whose major semi-diameters are equal to the earth’s polar 
and equatorial radii. Consequently, Eq. (2.3.6) may be expected to yield 
highly accurate results when used in an analysis concerned with a near- 
earth satellite idealized as a particle. 
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e = 0 



0 


Figure 2MA 


64 



2.9 


2.9 Differentiation with respect to a vector 

It is sometimes convenient to employ vector operations on scalar 
functions in order to express gravitational forces and moments. These 
operations may be interpreted as ordinary and/or partial differentiations of 
scalar functions with respect to vector variables. In some cases the vector 
in question is a position vector, and the ordinary derivative is then a 
spatial gradient . When the differentiation of a scalar function produces a 
force, the scalar function is called a force function . In other cases dif- 
ferentiations with respect to bther vectors are employed, as in the representa- 
tion of gravitational moments in terms of derivatives with respect to unit 
vectors. 

To unify the presentation of this subject in Sec. 2.10-2.20, the math- 
ematical tools there employed are first discussed briefly in this Section. 

If a scalar quantity F depends on a vector y, then it is useful to define 
a vector denoted by as follows: Introduce an arbitrary set of mutually 

perpendicular unit vectors a^, a^, and a^; let v^ 4 Y * ; regard 

F as a function of v^^, V 2 and v^; and set 


V 

-V 


- A . 3F ^ 9F 

£1 + -2 ^3 -3 


( 1 ) 


The vector F constructed according to Eq. (1) is invariant with 
respect to the choice of vector basis a^, a^, a^. The operation denoted by 
may be termed differentiation with respect to vector y. 

Similarly, if a vector quantity F depends on a vector v, and for some 
arbitrary vector basis a^, a^ one defines F^ A F • a^ and = y • a^ 
(i=l,2,3), then a scalar and a dyadic may be defined as 
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( 2 ) 


V 

™.v 



9F2 913 

9V3 


and 

ly I A (?v h) il * (^v ^ 2 ) H2*{\ ^ 3 ) i3 <« 

The quantities defined by Eqs« (2) and (3) do not depend upon the choice 
of vector basis a^^ a 2 ^ a^® 

In Eqs« (2) and (3) a vector differentiation concept is introduced 
without identification of a particular reference frame in which the differ- 
entiation occurs; this omission is possible only with the stipulation that 
the orientations of the unit vectors do not depend on v. 

Further useful quantities are at times obtained by cascading some of 
these definitions® For example^ when the operation defined in Eq. (2) is 

performed on the vector defined in Eq® (1) , the result is a scalar denoted by 
2 

¥t 

-v 

F A V • / V f\ (4) 

-t = \ -i / 


When y is the position vector £ of a generic point relative to some 

other point j and when the vector £ can be inferred from the context, the 

2 

subscript y = £ is often omitted; the quantities VF, V®F, VF and V F are 
then called respectively the gradient of F, the divergence of F, the gradient 
of F, and the Laplacian of F® 

As a consequence of the definitions given in Eqs® (l)-(4) and of 

various theorems of the differential calculus of functions of one or more 

2 

scalar variables, the quantities 7^ F, and V^F satisfy many relation- 

ships having counterparts in this calculus* For example, if v = w -I- c, 

^ ^ 
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where c is independent of and if F(v) and G(w) are functions such that 
G(w) = F(w+c)g then 

V F - V G (5) 

-V 


Or, suppose that v(gj^ , * , . ,g^) denotes a vector function of the scalar inde- 
pendent variables suppose, further, that F(v) and G(g^,..*,g^) are 

functions such that G(g^, g^) = F[v(g^, * , .,g^)] ® Then 



V F 

-V 



Cr*”!, e . « , n) 


( 6 ) 


Three differentiation formulas Involving a unit dyadic U and a unit 
vector u having the same direction as a vector v will prove useful in the 
sequel. These are 

V v - U (7) 

-V - - 

V V - u (8) 

-y 

and 


where v A y • u, so that v = v u« 

In dealing with gravitational moments, it is at times convenient to 
invoke the idea of partial differentiation with respect to a vector, in 
the following sense: suppose that a^, ^ 2 ^ a^, v, y, u and F have the same 

meanings as heretofore; let u^ A u * a^ (i = 1,2,3); and let y(u,v) and 
G(u,v) denote respectively a vector function and a scalar function of 


As a notational convenience, a scalar quantity such as F which depends upon 
the vector v is designated F(v); the functional representation of this 
quantity (as required by Eqs. (l)-(4)) is written 

F(vj,v^,vp, where ^ X * A Y ® (i=l,2,3) for arbitrary 

vector bases a^^a^ja^ and 2 x^ 22 ^ 23 ^ 
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Uj^,U 2 ,U 2 » and v, choosing G(u,v) such that G(u,v) ^FlvCujv)]. Then, defining 
9G/du as 


M -4 M X M 4. M 
3u —1 3 u 2 “2 3uj -3 


one can write 


V 

-V 



V u 

V - 



u 


( 10 ) 


( 11 ) 


or, in view of Eq, (9), 


V F . i M + /M _ i 
-V ~ V 3u \ 3v V 3u 



It follows from this equation and from v “ vu that 

V X V F = u X 
- -V - 3u 


( 12 ) 


(13) 


Derivations ; In characterizing a quantity as a function of a vector, 
one implies that the quantity is represented in terms of the scalar compo- 
nents of the vector by the same functional relationship, regardless of vector 
basis. Thus for example if the scalar quantity F is a function of the vec- 
tor V, and ^ V • a^^ and ^ v • aj^(i=l,2,3) for arbitrary orthogonal 
vector bases apa 2 ,a^ and bhen in functional notation 

F(Vi,V 2,V2) - F(v^,v^,vp (14) 

Using this relationship, one can prove the invariance of V^F, • F, 

2 ” ” 

V^F and V^F with respect to the choice of vector basis by establishing 

the equality of alternative representations of each of these quantities in 

terms of the arbitrary orthogonal vector bases a^,a 2 ,a 2 and apS^f end a^* 

This objective can be accompllshed'as follows; 
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Letting A (i»j“l»2,3), one can write 

Si ' sj <=« 

^ - 'j <=lj 


’i ■ '’J 'ji 

from the last of which It follows that 


9vj 31 

Consequently , using familiar differentiation theorems, one can write 


r. 9v! 3v. 3v* 

^ (14) ^ ^ (17) J 


3v! ^Ij 


^ a - C a ' C 

9v. -1 9v] ^Ij -k ^ik 


^ (18,15) ^ 


(1.2.15) 

Hence, 

9F 9F , . 9F , 

9v. -i “ 9v' 2k ^jk " 9v’ % 

(19,20) 3 1 


which establishes the Invariance of 7 F. For 7 * F, one has 

-V -V -* 


?. 9F' 9F' 9y/ 

— Si j r s 4 k o 

r 9v. ^ij 9# 9v. ^ij 

^ ( 16 ) ^ r ^ ^ 


(17) ^'^k (1.2.15)^ 
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and the proofs for F and • VF proceed similarly. 
When F in Eq. (3) is replaced with v, one has 

Now, 

9v 

^ 'V. = a, = 6. . a. = a. 

“Y ^ (1) “J 

Hence 


V V = a.a. = U 
-V - -i-i 

in agreement with Eq. (7). 

I 

If V = V u, so that u has the same direction as v, then 

V = (v • v) \ 


and 


V 

— V 


V 


= V (v*v)^ 
(24) " 



= 4 v“^ (v+v) 
(23,24) 


-1 

V V 


u 



( 21 ) 


( 22 ) 


(23) 


(24) 


which proves the validity of Eq. (8) . 

As for Eq. (9) , note that 

V u = V (v”^v) = ( V v”^\v + V ^ V V 
-V - “V - \-v y- -V - 

- (-V ^ 7 V I V + v~^ V V 

\ -Y J - -Y - 

= {-v^u)vu+v^U 
(7,8)' - 

= V ^ (U - uu) 

Finally, in connection with Eq, (11), note that both and v can be 
expressed as functions of Vj^,V 2 >V 2 (where v^^ A v * a-^)^ 
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(vj+v^+v^)-^,v.(v^+,|+v^) 


2 ^ _2 ^ 2\ls 


It follows that 


so that 


M » M ^ + M 3v 

9v^ "’ 9Uj 9vj^ 9v 


u w M 

-V * 9vj^ -1 ^2) 9vj^ -i 


( 25 ) 


9G „ . 3G „ 

" T— V u. + -5— V V 

( 1 ) ^ 



and this Is equivalent to Eq. (11) since u Is a syiometrlc dyadic, as Is 
apparent from Eq, (9). 

Example ; In Fig. 2.9.1, r,X and 3 are spherical coordinates of a point 
P, and hj^»h 2 »h 2 are unit vectors pointing respectively In the directions in 
which P moves when r,X,3 are made to vary, one at a time. The position vector 
£ of P relative to 0 may be expressed as 

£ (r,X,3) = r(cXc3aj^ + sXc3a2 + sBa^) (26) 

If F(£) denotes a function of £, and if G(r,X,3) is defined as 


G(r,X,3) A F[£(r,X,3)] 


then Eq. (6) can be used to express the gradient of F in the following 
frequently convenient form: 




1 9G , . 1 9G , 

rHF ^ ^2 + r ^ ^3 


(27) 
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To this end, one may begin by evaluating the partial derivatives of £ 
with respect to r,X and gs 


3£ 

•s— * cXcga- + sXcga- + sga, =■ b, 
(26) 

(28) 

_ . r(-sXcBa- + cXcga-) = rcgb„ 

(26) 

(29) 

3p 

^ ax r(-cXsga, - sXsga, + c3a_) == rb, 
93 ^26) 

(30) 

Next, from Eq, (6), and with V written in place of V^, 


ar 

p VF • -5=^ - VF • b, 

(28) - 

(31) 

1 i 

(32) 


(33) 

Now, since + ^ 2^2 + isi a vinit dyadic, 


VF - 7F • (bj^bj^ + b2b2 + b^b^) 


“ 2^ * -1-1 Z®" * -2-2 ^ * -3-3 

(34) 

Solving Eqs, (31)-(33) respectively for VF • bj^, VF • b 2 » 

and 7F * b^, and 


substituting into Eq, (34), one arrives at Eq. (27). 
The Laplacian of F may be expressed as 


V^F 


^2[3rf dv) 


, 2n 3 G , Q 9 / n 9G \ 1 


( 35 ) 


by using Eq, (4) in conjunction with Eq. (27), as follows: 


73 



2.9 


V F = V • VF = 
(4) (27) 


/ 3 . . 1 3 . 1 3 V \ 

\ 3r -1 rcB 3X -2 r 3& -3/' 


VF 


- Si • isl (ZW + 52 • -sf + 7 ^3 

Evaluating the partial derivative appearing in Eq. (36) but displaying only 
those terms that will not be eliminated by subsequent dot-^ltlpllcatlons, one 

has 2 

^ (VF) - ^ b + ... 

(27) ^ 

J. (VF) + J:_3!ib +iM!k + 

3X 3r 3X ^ rcB g^2 -2 * r 36 3X ^ *** 

3G Q, ^ 1 3^6 . 1 3G Q- ^ 

“ 3r *^^-2 r^E ^ -2 ■" r H *** 

~ 3r 36 r gg2 -3 *•* 

3G . .1 3^G . . 

” 3r ~3 ■*■ r ^^2 -3 


and substitution into Eq. (36) then yields 


„2_ 3^G ^ 2 3G ^ 1 

V F => — s- + — -r^ + 


3^G s6 3G ^ 1 3^G 


3r2 ■ - 9- ■ 3X^ " ^^06 ^ 


r^ 36^ 


which is equivalent to Eq. (35). 
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2*10 Force function for two particles 

The gravitational force F exerted on a particle P of mass m by a particle 
P of mass m (see Sec, 2,1) can be expressed as 

F - V V A W (1) 

- “E “ - 

where p is the position vector of P relative to P and V is given by 

V = top"^ + G (2) 

with p defined as 

V k{tf 

and C an arbitrary constant, while G is the universal gravitational constant. 

A scalar function of a vector variable is called a "force function" if 
the derivative of the function with respect to the variable is equal to a 
force. Thus V is a force function associated with the gravitational inter- 
action of two particles. 

The Laplacian of V (see Sec. 2.9) is zeros 

V^V = 0 (4) 

Eq. (4) is known as Laplace *s equation . This equation arises in many prob- 
lems of physics, such as electrostatics, magnetostatics, hydrodynamics, 
and heat flow, as well as in gravitation, and it is one of the classical 
equations of partial differential equation theory. Any solution of Laplace Vs 
equation is called a spherical harmonic, and in the context of gravitational 
problems a solution is classically called a gravitational potential . For the 
purpose of characterizing the interaction force F between two particles, 
however, only that special spherical harmonic given by Eq. (2) provides the 
force function appropriate for substitution into Eq. (1). 
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Derivations ; Differentiation of V with respect to p gives 

W = GimnVp = -Gtnmp u (5) 

(2) ■ (2.9.8) 

where u is the unit vector in the direction of p, that is, 

u “ p ^p (6) 

Consequently, 


VV = -Gmmp~\ = F (7) 

- ( 2 . 1 . 1 ) " 

■Ijhe Laplaclan of V is 

■ V^V = V • W =» -G5m7 • / p"^p ) 

I ' (2.9.4) “ " C7) ■ ' 

- -Gmm /sp ^u»p - 3p = 0 (8) 

(2.9.8, ' “■ ' 

2.9.2) 

Example ; A particle P of mass m and two particles Pj^ and P 2 , each of 
mass m, are situated as shown in Fig. 2.10.1, where r,X and 3 ate spherical 
coordinates of P and hj^»b 2 ,b 2 are unit vectors pointing in the directions in 
which P moves when r, X and 3 are made to increase one at a time. A force 
function V for the resultant gravitational force F acting on P is to be con- 
structed, and this function is to be used to express F in terms of components 
respectively parallel to bj^,b 2 ,b 2 . 

The forces Fj^ and F 2 exerted on P respectively by P^^ and P 2 can be 

expressed as 


2l 


V V 
-£i 


i 


(i-1,2) 


( 9 ) 


where p^ is the position vector of P relative to P^ and where Vj^(£^) ^nd 

V 2 (£ 2 ) given by 

V, » Gramp"^ + C. (i=l,2) 
i (2) i i 
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with equal to the magnitude of Hence 


F 


F 

-1 



+ V V 

“E2 


2 


( 10 ) 


If p is now introduced as the position vector of P relative to 0, then 

= p + c^ (i=l,2) where is independent of p^; if then W^(p) is defined 
as 


then 


Wj^(p) A V^(p + c^) (i=l,2) 


Vp = V (i=l,2) 

^ (2.9.3) E 


Hence 


F = V W. + ? W„ = V (W, + W-) = V V 

- (10,12) "E 1 -P 2 p 1 2^ -p 


( 11 ) 

( 12 ) 

(13) 


where V Is def ined as 


V A W, + W, = V, + V, = Gnm (p7^ + p"^) + C 

“ ^ ^ ( 11 ) ^ ^ ( 12 ) ^ 2 


(14) 


V is the desired force function. 

To express F in terms of components parallel to that 

Pi * ^r^ + - 2rLsin3^^ (15) 

and 

P 2 + 2rLsinB ^ ^ (16) 

and define a function W of r,X,3 as 

W A Gnm |^r^+L^-2rLslnB + ^r^+L^+2rLsing ^ j + C (17) 

Then 

V = W (18) 

(14-17) 
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and Eq, (27) of Sec. 29 provides 


, „ _ 3w , . 1 aw , . 1 . 

- 37 lx -2 7 H ^3 

= — Qnm I ^r^+L^-2rLsinBj ^r-LsinSj 


+ ^ r^+L^+2rLsing^ ^ r+LsinB jj 

f / 2 2 

Gram I (r^+L -2rLsin3j Lcos3 


-( 


r^+L^+2rLsin3 j 


-3/2 


Lcos3 


^3 


<19) 


Substituting from Eq. (19) into Bq. (13) one then arrives at the desired 
expression for F. 
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2.11 Force function for a body and a particle 

The resultant gravitational force F exerted by a particle P of mass 
m on the particles of a (not necessarily rigid) body B (see Sec, 2,2) can be 
expressed as 

F « (1) 

A _ 

where R is the position vector of the mass center B of B relative to P 
and V is a force function given by 
N . 

V = ® C <2) 

1-1 ^ ^ 

with 

Pi=(e^)^ (1=1,,,., N) (3) 

and 

PiA5 + r^ (1-1,..., N) (4) 

Here B Is presumed to consist of particles F^,..., F^ of masses m^,..., 

mjj, Pj^ is the position vector of F^ relative to F, assumed nonzero, r^ is the 

A 

position vector of F^ relative to B, C Is an arbitrary constant, and G Is 
the universal gravitational constant. If B Is a continuous distribution of 
matter not including the point occupied by P, then V Is given by 


with 


and 


V 



pdT + C 



(5) 

( 6 ) 


p = R + r 

where p is the mass density of B at a generic point P of B, p is the posi- 
tion vector of P relative to P, r is the position vector of P relative to 
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* 

B , and dx is the length, area, or volume of a differential element of the 
figure (curve, surface, or solid) occupied by B. 

The force functions in Eqs. (2) and (5) both satisfy Laplace’s equation: 


V^V = 0 

Derivations; Differentiation of V with respect to R gives 
N _ 

V V = -Gm X) m p"'^ V p 
5 (2) i=l ^ ^ i 


(7) 

( 8 ) 


Now, 


Vi = 

- ^ (3,4) - 






21 




(Rfr^) 


(3,4) 


Pf Ei * \ ^^+£l^ = P^^Pi * (U + 0) = 


( 9 ) 


Hence 


N 

V V = -Gm 2^ “iPl^Ei = I 
- (8,9) 1-1 (2.2.1) 


A parallel derivation shows that Eq* (1) remains valid when Eq. (2) 
is replaced by Eq, (5)« 

Eq, (7) follows directly from Eq. (2,10.8), since the order of the 
operations of evaluating the gradient and integrating over B is interchange- 
able. 

Example s A uniform, thin rod B of length 2L and mass m is subjected to 

the gravitational attraction of a particle P of mass m. The value of a force 

function V associated with the resultant gravitational force exerted by P on 

B is to be determined in terms of the cartesian coordinates x,y of P for a 

* 

coordinate system with origin at the mass center B of B, as shown in Fig. 
2.11.1. The constant C in Eq, (3) is to be chosen such that V approaches 
zero when y=L and L/x approaches zero. 
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If the rod Is regarded as matter distributed along 
segments then the distance p from f to a generic point P 
expressed as 

f 2 2 ^ 

p = (x + Ti ) 

where Ti varies from y~L to y+h. The mass density p of B 
m/(2L), and a differential element of B has a length dn* 


V 


(3) 



A 


dn + c 


and, when y«L, 


Gmm I yHL + + 


2L fi 

y-L +>/x + (y-L)“ J 


+ G 


itim V 


L 


0 


X 

which equals zero 


SS C 


by hypothesis. 


a straight line 
of the rod can be 

at F Is equal to 
Hence, 
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2»12 Force function for a small body and a particle 


When a body B is subjected to the gravitational attraction of a particle 

ma ^ 

P removed so far from the mass center B of B that the largest distance from 
* 

B to any point P of B is suhstantially smaller than the distance R between 
— * 

P and B , a useful- form of the expressions for the force function V given in 
Eq, (2elle2) and (2,11,5) can be found as follows: Replace p^ and p with 

|^(R + with |^(R + r)^ ^ , respectively j and then expand the integrand 

in ascending powers of ]r |/R or |r|/R to obtain 


Gmm 

R 


(‘ • I. ■“’) 


+ c 


( 1 ) 


where is a collection of terms of the i^^ degree in |r^|/R or [r|/R, m 

and m are the masses of B and P, C is an arbitrary constant, and G is the 

.( 2 ) 


th 


universal gravitational constant. In particular, v 

„( 2 ) _ 1 , 


is given by 


2mR^ L 


tr(I) - 31 


11 


( 2 ) 


where tr (I) is the trace of the inertia dyadic of B for B , and I is the 

* 

moment of inertia of B about the line connecting P and B , so that (see 


Fig. 2.3.1) 


111 A • 1 ‘ 3i 


(3) 


Because depends upon the orientation of a^ relative to B, Eq. (2) 
is sometimes less convenient than an alternative form involving central moments 

f ? I 

and products of inertia of B for an arbitrary vector basis here 


represented by 

I 






(4) 


The required expression is 
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[41 




(^12‘^iA 2 ■'■ ’■13*^11*^13 ■*■ ’23*’l2*^13/ 


*=13 ^ Si • 

f t I 

When B Is a rigid body and ^3 fixed in B, the scalars 

(j ,k=l5 2,3) in Eq* (5) become constants, whereas in Eq* (2) remains a 
variable. 

For the special case in which the body-fixed unit vectors parallel 
principal axes of B for B , past conventions (see Sec. 2.3) suggest the nota- 


.(2) _ _1_ 


j^Il ^ 1 - 3 C^^) + I2 (i-3Cj^2) ^3 { ^“^^13) 


where and (j=l,2,3) are defined respectively by Eqs. (2,3.8) and (2,3.9), 

Eqs, (5) and (6) adopt particularly useful forms when the Indicated 
direction cosines are written In terms of the spherical polar coordinates 
shown In Fig, 2.12.1. Eqs. (5) and (6) then become respectively, 

^ - 2I33) ( 3 <ll>>Vl)- 3 (lii - cos2scos2X 

-6 ®o*^Psln2X + l^^sin28cosX + I23 sin28elnx)j (7) 

and 

^ [(’1 ’2 - -1) (8) 


+ I2 - 2l2M3sln 3 




cos pcos2X 


It is often convenient to introduce some measure of the dimensions of 


B into v^ in order to obtain an expression involving dimensionless 
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Figure 2.12.1 
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constants and dimensionless variables. If the symbol Rg Is assigned to the 
selected normalizing dimension, and the symbols 


J A 


-1 


and 


e A - 


2mR 


2 »1 + ^2 - 


2 

B 


(9) 


( 10 ) 


are introduced, then Eq. (8) becomes 

2 


(sslnVl) + I (^) cosWs2X (11) 

If F is the approximation to F defined in Eq. (2.3.6), then F can be 
expressed as 

( 12 ) 


I " 


where V is given by 


^ &nm 
R 


( 1 + v^^^) + G 


(13) 


and R is the position vector of B relative to P. V' 

In other words, V is a force function suitable for dealing with the 
gravitational interaction of a body and a remote particle when one uses F to 
approximate F, 

As previously, the force function V satisfies Laplace’s equation; more- 
over, V also satisfies Laplace’s equation. 

1*5 


r 2P 

Derivations t Replacing p in Eq. (2,11,5) with j (R + r) | provides 


V ^ dm J ^ 


+ 2R • r + r^l p dx + C 


(14) 


Introducing q as q ^ r/R, noting that R/R = and applying the binomial 

2 

series, convergent for |2R«r + r | < 1, find 
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j ^X+2aj* q+t^'^ pdT + C 

I I {^ll) + •••]p dT + C 

|i jy*p dx - ^ .y r p dx - ^/[E^-3(a^.r)2] pdx 


+ ... + C 


p dx = m 


J" r p dx = 0 

y ^ P dx - 3 y|^^-^aj^»r^^p dx -2 ^ 


Furthermore, 


Hence, 


y"p- (si'£ ) ] p dT = ij 

2 J“ P dT “ tr(l) 

/ "" ( 18 , 19 , 


2 / r p dx =s tr(l) 


20 ) 


V = Gmml 1 + Itr(l) - 31, ,1+ ... ■+ C 
(15,16,17,21) R 1 2mR L " J 

In agreement with Eqs. (1) and (2). 

Eq. (5) follows directly from Eq. (2) with the substitution of 
Eq. (3) for the (summation convention) expansion 

i (I) 
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and the definition of . With these relationships, Eq. (2) becomes 


(2,3,4) 2mR 


^ (^11 - ’5l • 5l) 


1 / » I f I \ 

“ 712 (\i " ^ ^li ^ij ‘^Ij ) 


“ [^11 (^-3Cn) + (1-3^12)-^ ^3 (i-3<3) 


( ^ 12 ^ 11^12 ^ 13 *^ 11^13 ■*■ ’• 21 ^ 12^11 ■*■ ^ 23 ^ 12^13 




32^13^12)] 


With the symmetry relations = 1^^, Eq, (23) confirms Eq, (5). Eq. (6) is 
merely the special case of Eq. (5) with vanishing products of inertia and new 
notation, and Eq. (7) and (8) are the special cases of Eqs. (5) and (6) in 
which the relationship 

I » 1 

a^ » cosScosXbj^ + cos8slnXb2 + sin3b2 (24) 

from Fig. 2.12.1 has been used to obtain for the required direction cosines 
the expressions 

» » 

Cii = cos3cosX ; ® cos3sinX ; “ sin3 (25) 


The combination of Eqs, (5) and (25) produces 


^(2) ^ 


[^n( 


l“3cos^8cos^X j + I22 (l- 


2 2 \ 
3cos 3sln \f 


’ / 2 \ / ’ 2 ’ 

+ 1^2 ( l-3sln 3y -6 I Ij^2^‘^os 3cosXslnX + I^^cos3cosXsln3 

+ I 22 sin3cos3slnXj^ 

- -JL , [l-3cbs^3 (Jg + % cos2X)j + 1*2 |^l-3cos^3 ih - h cos2X)j 


89 





2.12 


6 


^ luu _ _ niL , .. , 

tr(I) = -j- j I,, = a, » I ® = 1 j- sin ^ 


11 =1 * 

after which substitution into Eqs. (2) and (13) gives 

.2 


~ Gmm 

Y SS . \ 

R 


1 - 


24R 


( 2-3sin^tJ;) 


+ C 


( 31 ) 


In order to obtain this result by means of Eq. (11) rather than Eq. (2), 
one must recognize by comparison of Figs. 2.12.1 and 2.12.2 i:hat in this 
example 0 = 0 and X - Substituting this result together with the mass 

properties 


T T “L 

h “ ^3 “ 12 


into Eq. (11) provides 




(32) 


(33) 


where 


j = 

(9) 2mRg ' 


0 + 


T 2 o T 2 

ml 2mL 


12 


12 


) 


24R, 


(34) 


and 


(10) 2mR 


2mR„ ' ' 

B 


(35) 


Although it is customary to make some specific choice for R„ (such as 

15 

R^ = L/2), in order to obtain numerical values for J and e, it is apparent 
that in combining Eqs. (33 - 35) one can cancel Rg fro,m v' so the choice 


of Rg is irrelevant to v 


( 2 ) 


,( 2 ) 


(33-35) R" 
l2 

“ —2 
48R-^ 


The result then becomes 


1 J 1 OI 


1 + 3 


^l-2sin^i|^y| = ^ ^4-6sin^iJ»^ 


(36) 
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Substitution of Eq. (36) into Eq. (13) confirms Eq« (31), as assured 
by the equivalence of Eqs. (2) and (11). For the orientation shown in Fig. 
2.12.2 the calculation from Eq« (2) is more direct, but Eq. (11) offers the 
advantage of explicit separation of mass properties and general kinenatlc 
variables, all in terms of dimensionless quantities. 

To verify that Eq. (12) is indeed satisfied, one may proceed as follows s 




- (31) R 


Gmm 

2 


V„R fl + f2-3stn^i|()l 

- L 24 R ^ ' J 


% [- ^ If - ^ V] 


Now, 


Hence 


R = R ^ - sln'|'b2j 


Vg^R = Vj^R ^cosil^bj^ - slnij^b2^ 

+ R l^-sinif) bj^ - coBijj V^\p ^2^ 

or, In view of Eqs. (2^10,6) and (2.10.7). 

U = ^cosijibj - sln?J;b2j - R V^ijj ^sinijib^ + cosi(ib2j 

“ 2i2i - R V 22 

Dot-multiplication with ^2 thus gives 

«2 ■ -R If 

and using Vj^R “ a^^ and one obtains 

2 2 
+ ( 2-3sin^ii^ ) a. - — , sin^ij) 




(37) 


(38) 


94 



2a3 


2^13 Force Function in terms of spherical harmonics 

When a body B of mass ^k is subjected to the gravitational attraction 
of a particle P of mass m removed so far from the mass center 8 of B that 

k 

the largest distance from B to any point P of B is substantially smaller 

» k 

than the distance R between P and B ^ the force function V given in Eqs, 
(2*11*2) and (2.11.5) can as an alternative to Eq. (2.12,1) be represented 
in series form as 


V 


Gmm 

R 




+ C 


( 1 ) 


where G is the universal gravitational constant, C is an arbitrary constant, 
R„ is an arbitrary distance such that R_ < R, and 


CiA 



( 2 ) 


with r representing the magnitude of the vector from B* to a generic point P, p 
representing the mass density function of and with P^(cosa) representing 
the Legendre polynomial with argument the cosine of the angle a defined by 
Fig. 2.13.1. In this figure , R,Aj8 are the spherical coordinates of P, 
r,y,Y the sperical coordinates of P^ and a is the angle subtended by the 

^ k 

two lines joining P and P to B . The series in Eq. (1) converges for (r/R)<l. 
By definition, the Legendre polynomials are given by PQ(y)Al and 


?i(y) 4 


1*3» 


(2i-l) 

if 



1(1-1) 1-2 . 1(1~1) (i-2)(i-3) i-A , 

(2i-l)«2 y (2i-l)(2i-3)«2*4 ^ 


. 


(1=1,..., oo) 


(3) 


with the bracketed series terminating with the y° term if i is even and 
with the y^ term if i is odd. Thus 

PQ(y) = 1 (^a) 

P2^(y) “ y (^^^) 
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Fiflurs 2.13.1 
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P2<y) = % (3y -1) (4c) 

P 3 (y) - h (5y^-3y) (4d) 

Eq. (2) for (1=2,..., <») can be written in such a way as to separate 

the properties of B from the location of P relative to B by introducing the 
associated Legendre function of the first kind, P^(y), defined by 


2^i!(l-j)! 


-j _ (fej)(lrj-l) 

(21-1) 2 y 


+ ^i- jJ!ja r i ri)a- j - 2),, (i-j-3) i-j-4 1 0 ^ . 

(21-l)(2i-3)»2«4 y + ...J U,3 u,i,..., ^ 

with the bracketed series terminating with the term if i-j is even, and 
1 

with the y term if i-j is odd. For specific indices of particular interest, 
pj(y) = Pi(y) (6a) 

and 

P 2 (y) = (l-y^)^ ; P 2 (y) “ 3y(l-y^)^ ; “ 3(l-y^) (6b) 

In terms of associated Legendre functions, 

= 53 (sin3)(G^, cosjX + S^.slnjX) (7) 


where 


^i0 4 


r A 2 (1-1) 
^Ij - m Ti+jT 




( 1 = 2 ,..., ~) 


(siny) cosjp pdT 


/i=2,...,oo\ 
\j=l,...,<» / 


S A ^ 

^ij ^ m (i+j) 




(sin3) slnjp pdx 


/ i“2, \ 

\j=l,...,«>/ 


so that the constants and (for all Indices) establish the properties 
of B without reference to the location of P. 
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When Eqs* (1) and (p are written in combination as 


V = 


Gmxn 


1 -h 


oo i / R \ i 

Z £q(r^) 4 (singHG^j cosjX + sinjX)] + c 


( 11 ) 


the force function ¥ is said to be in spherical harmonic form. Terms such 

as (sinB) cosjX and P^ (sing) sinjX are classically called tesseral 
th th 

harmonics of i degree and j order; tesseral harmonics of order zero are 

th 

called zonal harmonics . A linear combination of tesseral harmonics of i 

degree is called a surface spherical harmonic of i^^ degree, and the product 

of R with an 1^^ degree surface spherical harmonic is called a solid 

th 

spherical harmonic of i^ degree. Thus Eq. (11) is a representation of V in 
terms of solid spherical harmonics. The International Astronomical Union 
has established Eq. (11) as a standard representation of the force function 
of the earth, choosing the mean equatorial radius of the earth for R^. This 
series converges only outside of a sphere containing all of the mass of B, 
since convergence of the series in Eq. (11) requires r < R for all r; due to 
the oblateness of the earth, the series does not converge for points in a 
region of space above the surface of the earth ranging in thickness from zero 
at the equator to approximately 20 km above the poles. 

In comparing the alternative series expansions for V represented by 
Eqs. (2.12.1), (1), and (11), it is helpful to recognize that 
1 


^2 2 
2mR" 


tr(I) - 31 


11 


( 12 ) 


where I is the inertia dyadic of B for B , and is the moment of inertia 

jfe mm 

of B about the line joining B and P, The equivalent contribution to Eq. (11) 
is given by 
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^2 = 1*2 cosjX + S 2 J sinjX) 

= h ( 3 sin^g"-l) C2Q + ^ sin 2 g (C22COSA + S2j^sinX) 


+ 3 cos 3 (C 22 GOS 2 X *f S22^i^2X) 


( 13 ) 


Physical interpretation can be provided for the constants and 
S 2 J for j = 0,1,2 by aligning the unit vectors bpb^, and b^ in Fig. 2,13.1 
with principal axes of B for B ; then 

1 


^20 “ , „2 ^^1 ’•2 " ^^ 3 ^ 

znUx— 

x> 


^21 ®21 “ ^22 “ ® 


(14a) 

(14b) 


and 


^22 


4 inR ! 


2 »1 - ^2> 


(14c) 


where 1 ^, and are the principal moments of inertia of B for B 


More generally, when and b^ are arbitrary orthogonal unit 

vectors, then 


20 “ ^^11 ^22 " 

2 l ^ 3 )/^ 2 mRg ) 

( 15 a ) 

21 “ "^ 13 ^ ("* b ) 

» ^22 “ “^^11 "■ ^ 22 ^^(^“^ b ) 

( 15 b ) 

'21 “ "^ 23 / ( “^) 


( 15 c ) 


For an axisymmetric rigid body B with its symmetry axis aligned with 
b^ in Fig. 2.13.1, the constants (i, j“0, « . . ,<») and (i«2, . . . 
j=l,...,<») in Eq. (11) are zero, and the remaining cons^aiitk C^q are tradi- 
tionally relabeled such that 

Cio 4-^1 (16) 
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Then V takes the form 

m (T 

Gvm 


V = 


R 


00 

1=2 ^ 


( 3 ‘ 


P^(sin3) 


(17) 


in which only zonal harmonics appear. 

In the notation of Sec. 2,12 (see Eq, 2,12.13), V may be approximated 
in some cases by 

Gw 


R 




+ C 


Laplace's equation is satisfied by both V and V, so that 
V^V = 0 


and 


2~ 

V V = 0 


Moreover, Laplace's equation is satisfied by 

>i+l 


i /I \ 1 r 

pJ (sing) AjCOsjX + sinjX 


(18) 


(19) 


( 20 ) 


( 21 ) 


for any constants and so that any truncation of the series expression 
for V in Eq. (11) consisting of linear combinations of terms given by Eq. 

(19) is itself an approximation of V which satisfies Laplace's equation. 

Derivations; Application of the law of cosines to Fig. 2.13.1 provides 


2 2 2 
p = R - 2rRcosa + r 

permitting the representation 


( 22 ) 


V = Gm 
(2.11.5) 




pdT + C 


( 22 ) 


Gm 

R 




pdx 
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The identity 


(l-2xy+y^ ^ y^P.(x) 

' i=0 ^ 

permits the force function to be expressed in terms of Legendre polynomials 
(see Eq. (3)) in the form 


. 00 
Gto 

R 




(cosa) pdx + C 


(23) 


whenever 


fl<i 


for all r, so that the distance r from B to P is substantially less than 
R for all P in B* The first two terms in Eq* (23) simplify, since 


and 


so that 


P^ (cosa) = 1 

^ (4a) 


Pj^(cosa) « cosa 


/(t)° 


-/ 


Pq ( cosa) pdx = / pdx = m 


(24) 


(25) 


(26) 


and by virtue of the definition of mass center B 


/(t/ 




Pj^(cosa) pdx “ / rcosapdx ® 0 


Hence 


(23,26,27) 


Gmm 

R 


[l+if /(f) Pj(C05C.)pdT] 


+ c 


(27) 


(28) 


and, by introducing the arbitrary Rg with the identity 


See O.D. Kellogg, Foundations of Potential Theory . Ungar Publishing Co,, 
New York, 1929, p. 128. 
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substituted into Eq® (28)^ one finds Eqs® (1) and (2) conflTmed. 

Eq. (3) is a definition^ which leads directly to Eqs. (4) in special 
cases® Similarly the definition in Eq« (5) leads directly to Eqs. (6)* 

Proof of Eq» (7) requires the trigonometric identity 

cosa = slnPsiny + cos3cosYcos{X-y) (29) 

which follows from the law of cosines applied to the spherical triangle in 
Fig. 2.13.1 defined by points B*, P, Q, and Q. With Eq. (29) and the lemma* 


i 

P . (cosa) = P , (sin3) P . (siny) + 2 Y*. 

i i i 


(izii! 

(i+3 ) t 


P^(sinB) P^Csiny) cosj(X-y) 

(30) 


Eq. (2) becomes 




P. (sing) r / vi 

,30) 5T“ J 


+ f TSd)f Pl(«i“6)j^cosjX 

+ sinjX J P^(siny) sinjypdx 


/ (s)'i 


(siny) cosjypdT 


(31) 


confirming Eqs, (7) - (10). 

Eq. (11) is the direct combination of Eqs. (1) and (7). 
Eq. (12) follows from Eq. (2), since 


*W.E. Byerly, An Elementary Treatise on Fourier Series and Spherical, 
Cylindrical, and Ellipsoidal Harmonics , Ginn and Co., Boston, 1893, p. 211. 
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=2 y(i^f 

= \ f r^(3cos^a-l)pdT 

(4c) 2mRg J 

= ^ 2r^pdx - 3 ^ r^sin^apdx J 

niRg 


where tr(I) is the trace of the inertia dyadic of B for B and is the 

ifc mm 

momeAt of inertia of B about the line joining B and P, 

Eq. (13) is a special case of Eq. (2), with Eq, (4c) and a pair of 
trigonometric identities substituted, 

jThe physical interpretations in Eqs. (14) and (15) are confirmed by compar- 
ing Eq. (13) with Eqs. (2.12.8) and (2.12.7), noting that 


.(i) 


(2.12.1,1)' 


!W‘ 


i=2,.,.«) 


(32) 


and equating coefficients of independent trigonometric functions, such as 
(3sin^3-l) , cos^3cos2X, cos^3sin2A, sin23cosX, and sin23sinX. 

Eq. (17) is a special case of Eq. (11) in which the definition in 
Eq. (16) has been substituted and due to syiranetry all dependence on X has been 
eliminated. 

( 2 ) 

Eq. (18) follows from the substitution of from Eq. (32) into Eq, 


(2.12.13). 


Eq. (19) is identical to Eq, (2,11,7), since V as expressed in Sec, 

2,13 is equivalent to V as defined in Sec, 2,11, 

The statement in Eq, (20) that V also satisfies Laplace’s equation was 
asserted but not proven in Sec. 2,12. Comparison of Eqs, (18) and (13) reveals 


that 
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V = 


R 




R / S P^Csing) (C2^cosjX + S2jSinjX) 


(33) 




so that V may be written as a linear combination of terms such as defined 
by Eq, (21). Thus by superposition Eq. (20) is valid if satisfies Laplace’s 
equation. 

In order to prove that is a particular solution of Laplace's equation 
when written in terms of spherical coordinates, note that V^v =. 0 becomes, 
from Eq. (2.9.35), 


i_ F-i- /r2 

^2 [ 9R dR) 


^ 2„ a^v ^ Q a 

+ sec p s- + seep 

ax^ 


Hw)] = 


(34) 


To obtain a general solution by separation of variables, adopt the assumed 
solution 


V = X(R)Y(g)Z(X) 

2 

and multiply by R cosg/ (XYZ) to obtain 


cos p 


d /T,2dX\, cpsg d f_ _r,dY 
dR dR / 






d^Z 

dX^ 


(35) 


(36) 


Since the left side of Eq. (36) does not depend on X, neither does the right 


side, so that both sides equal the same constant, here called j . Hence 

~ + j^Z = 0 
dX^ 


(37) 


and 


cos B d 
X dR 




Eqi (37) has the solution 


Z = cosjX + B^sinjX 
(37) - ^ 


(38) 


(39) 


for arbitrary constants and 
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Eq. (38) can be divided by cos'^6 and rewritten as 


L.J. /R2dx\ 
X dR dR/ 


1 111 odY 

(3« 


The left side of Eq. (40) depends only on R, and the right side only on 3; 
hence both must be a constant, here designated 1(1+1). Thus 




i(l+l)X 




Eq* (41) has the solution 

X = A-R^ + 

(41) ^ ^ 

for arbitrary constants A2 and as may be confirmed by substitution. 

Eq. (42) adopts a more familiar form when written in terms of 


y A sin3 


so that 


( )' 4 { 


cos3j^ 

dy 


The resulting equation is 




+ i(l+l) - — 5- Y - 

L dV) J (42,44) 


(lV)^-2 y^+[i(i+D I t . 0 

dy I (l-y)J (45) 


Eqs. (42), (45), and (46) are known as associated Legendre equations 
of degree 1 and order J. If the symbol F^(sln3) Is defined to represent a 
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particular solution to Eq. (42), then Eq, (34) has a particular solution 
=s (sin3) fA^cosjA + Bj^sinjAl 


(35,39,43) ^ L JL J 

which implies a second particular solution 

i^l 1“ T 

^i " (r) p{(sinB) AjCosjX + BjSlnjX (47 

thus proving the contention that as defined by Eq, (21) satisfies 
Laplace’s equation, subject to the implied supposition th^t the function 
P^(sin3) satisfying the associated Legendre equations is the same as the 
P^(sinS) defined by Eq. (5). 

i 

In order to confirm the implied meaning of P*^(sin6), it is convenient 
to rewrite Eq. (46) in terms of W, defined by 


so that 


W A 


Y = 

(i-y ) w 

11 

li-v 1 37- 

A 

dy^ 

- 

dy 






j/2-l 


Substituting Eqs. (49) - (51) into Eq. (46) and multiplying by ( 1-y ) 
yields 


2 \ 2 


li 2\ d W 

( 1 -y ) -3 

dy 


2y(j+l) || + [i(i+D - j^-j] W = 0 


Eq. (52) is more meaningful when written in terms of L, as defined by 
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is a particular solution of Eq. (46), In order to confirm Eq. (5), it is now 
necessary first to establish Eq. (3) for Pj|^(y). To this end, let 

aA(y^-l)^ (57) 

SO that 

^ = 2iy ( y^-l) 

and 

( = 0 

Differentiate once with respect to y to obtain 

(l_y2 ) + 2y(i-i) || + 2ia = 0 

dy ^ 

and differentiate with respect to y i times more to establish 
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L 2\ / d^a\ 

^ ^ 'dy^ \<iy / 


2y 


d / d g 

\dy^ ; 


+ i(i+l) 


($)■• 


proving that d^a/dy^ satisfies Legendre's equation. Since Pj|^(y) has been 
identified hefe simply as a particular solution of Legendre's equation, one 
can define P^(y) as an arbitrary multiple of d^a/dy^, and specifically let 

\± 


Pi(y) A 


(57) 2-"i! dy 




(58) 


This result, known as Rodrigues* formula for the Legendre polynomials, is 
usually considered as a property of the particular solution satisfying 

Legendre *s equation, but a systematie evaluation of the derivative in Eq. (58) 
confirms the equivalence of this expression and the series developed in 
Eq. (3). Substitution of this series into Eq. (56) and repeated differentia- 
tion confirms the series expansion for P^(y) developed in Eq. (5). 

Example ; In conventional practice Eq, (11) is used in application to 
natural celestial bodies for which mass distribution properties are not known 
sufficiently well either to calculate and from Eqs. (8) - (10) or to 
use the series for V pre'sented in Sec. 2,12, The alternative is to deduce 
values for and S^^ (at least for the lower values of i and j) from 
empirical data on gravitational forces on other bodies. 

In the case of the earth, for example, measurements of gravitational 
forces at the earth’s surface and on artifical earth satellites have estab- 
lished the approximate values 
= -1.0827x10"^ 

% = % = 0 = i»2) 

when the normalizing quantity is given the value of the earth’s mean 
equatorial radius, here designated b. 
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If the earth is assumed to be a homogeneous spheroid, with polar radius 
G, the empirical data can be used to compute the ratio c/b. Since 

« I2 - m(b^+c^)/5 

and 

13 - 2mb^/5 


we have 


-1.0827x10 


-3 


(14a) 2mb 


1 r 2m 

K?-‘) 


24 ^ 2 ) _ ^ 


and 


£ 

b 


(1-5x1.0827x10“^) 


H 


0.9973 
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2,14 Force function for two small bodies 

^ 

When the distance R between the mass centers B and B of two (not 
necessarily rigid) bodies B and B exceeds the greatest distance in either 

•Si > ■ 

body from the mass center to any point of the body, so that the resultant 
of the system of gravitational forces exerted on B by B can be approximated 
by a force F defined as 

+ +£®) (1) 

where G is the universal gravitational constant, m and m are the masses of B 

- , ( 2 ) -( 2 ) 
and B, a^ is a unit vector directed from B to B, and f and f are given 

rs^ 

by Eqs. (2.4.2) and (2.4.3), then a force function V such that 


I where R = Ra^ , is given by 

C; Gm (2) . -<2) \ 

V = -g — (l + v '+v l + c 


( 2 ) 


(3) 


with 


and 




|\r(I) - 31^ J 

l^tr(f) - 


(4) 

(5) 


where tr(I) and tr (I) are the traces of the inertia dyadics of B for B and 
of B for B , while and are the moments of inertia of B and of B about 

it 

the line connecting B and B . 

Derivation : Differentiating Eq, (3) with respect to R and then proceed- 

ing as in the derivation of Eq. (2,12.13) , one arrives at 
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and use of Eqs. (2.4.2) and (2.4.3) then gives 

i 

V V- -2% (a, + +£«>) - F 

Example ; The Example in Section 2.4 deals with an approximation F to 
gravitational force exerted on a rectangular parallelepiped by an oblate 
spheroid. To construct a force function V that satisfies Eq. (2), note that 
the moments of Inertia and can be expressed as 

111 - !r 1 <si- V' 

ill - # ^ (51- ii>^ i^] 

while 

2 

. /t*\ root /o r-.^\ 

tr (I) = ~g- (3 - e ) 

and 

tr (T) = (3 - i^) 

V can then be formed by substituting into Eqs. (3) - (5); 



and it may be verified that differentiation of V with respect to R (using 
Eqs. (2.9.8) and (2.9.9) leads to F as given by Eqs. (2.4.14), (2.4.26), and 
(2.4.29). 
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2 * 15 Force functions for a centrobaric body 

Since for purposes of dealing with gravitational effects a centrobaric 
body (see Section 2.5) may be replaced with a particle situated at the mass 
center of such a body, force functions suitable for dealing with the gravita- 
tional interaction of a centrobaric body and a particle, a centrobaric body and 
any body whatsoever, or a centrobaric body and a remote body can be obtained 
directly from Eqs. (2.10.2), (2.11.2), (2.11.5) , (2,12.1) (2.13.1) or (2,13.11). 

Example ; The gravitational force F exerted by an uniform sphere S of 
mass M on a dumbell of mass 2m (see Fig. 2,15,1) can be expressed as 



with 

V - GMm 

( 2 . 12 , 2 ) 

where n is a unit vector directed as shown in Fig, 2,15.1. 


|^(R + + |^(R - Ln)^ j 




+ C 
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m 



Figure 2.1S.1 
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2,16 Force function for a body and a small body 

When the distance between the mass centers B and B of two bodies B and 

* 

B exceeds the greatest distance from B to any point of Bj and a force function 
V(g) for tl>e force exerted by B on a particle of unit mass at a point P situated 
as shown in Fig. 2.16.1 is available, then there exists a function V(R), where 
R Is the position vector of B relative to B , such that the resultant of all 
gravitational forces exerted by B on B can be approximated by a force F expressed 
as 

F = VV(R) (1) 

/ 

wherd the synbol V connotes differentiation with respect to R. The function 
V(R) is given by 


I', V(R) = mV(R) - f I: VW(R) + C (2) 

where m is the mass of B, I is the inertia dyadic of B for B , and C is an 
arbitrary constant. 

Derivations t The gravitational force dF exerted by B on a differential 
element of B at P can be expressed as 


dF = V VpdT 

( 2 . 11 . 1 ) “2 


(3) 


where p is the mass density of B at P while dx is the volume of the element, 
dpnsequently, the resultant F of all gravitational forces exerted by B on B 
is given by 


F = y* VpVpdT 


( 4 ) 


The double dot product in Eq. (2) is defined such that, for two dyads 
and v^Y 2 » (u^U 2 ) : (YiYy^ “ -1^ ^-2* -2^ • obeys the distributive 

law when applied to dyadics; 

(a^a2+b^b2+. •.): (A^A2+Bj^S2+. . • ) = ^±^ 2 " -1-2 ^ -1-2 * §1-2 


*4* bj|^b2 • ^ ^1“2 * —1^2 • • • 

+ • • • 
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Figure 2.16.1 
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To obtain the desired approximation F to F, one can expand V(£) in a 
Taylor series about g ® R, retaining only terms up to and including those of 
degree three in jr|, where r A g - R (see Fig. 2.16.1); differentiate with 
respect to g; and then carry out the integration indicated in Eq. (4), These 
tasks are facilitated by introducing a function W(r) ^ V(r + R), which permits 
one to write 

V V = V W (5) 

"2 (2.9.5) 


and 

F = /*V WpdT (6) 

" (4.5)*' 'S 


If mutually perpendicular unit vectors and ^ r • n^ 

(i=l,2,3), the Taylor series expaUsion of W about r = 0 can be written (using 
the summation convention) 


W = W(0) + r, 




1 

2 ! 




W 




3! 




W 


,ijk 


(7) 


where 


W ^ A 


9W 


,i - 3r, 


3^W 


r=0 


(i»l,2,3); A 


(i,j = 1,2,3) 


r»0 


and so forth* Differentiating with respect to r^, one obtains 

= 0 + 6..W . + 6,»r, + r.6,» ) W . . 

(7 8) 21 \ iil 3 i jJl/ ,ij 


+ ’^i'^jJl’^k + ’^i’^j-kJl J ",ijk 




+ * , , 




Hence, 


( 8 ) 


V w = /V W\ + r • /V V W\ + i r r: /V V V w\ + 

“E (2.9.1) TH ; ^ - \-l-l I r. / n 

r=0 r=0 r=0 

Substituting into Eq. (6) only the terms here displayed explicitly, and 

noting that 
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/pdT ® I& ^ y* rpdx » 0, ^y* crpdT ■ ^ tr (1) U - I 
where U is the unit dyadic, one can now write 

? - m(j b) + T [t “ <P H - ll = ( XrWr’') 

' - 'r=»0 J \ 

or, after using Eq. (5), 

? ■ (Yj, + 1 [i “(P H - 1] ! ZjYI * °) 


g"R 


where C is an arbitrary constant. Now, differentiating V(g)with respect to g 
and then setting g equal to R Is precisely the same as differentiating V(R) 
with respect to R. Hence, If V connotes differentiation with respect to R, 
then 

F = V j mV(R) + ^ tr (I)U - I j s 2W<R) + C j 

Now 

U; VVV(R) » V»VV(R) » V^V(R) » 0 

(2.11.7) 

Thus Eq. (1) follows immediately if Eq. (2) Is used to form V(R). 

Example ; When V(R) is available in the form of an explicit function 
V (R,X,3) of the spherical coordinates R,X,3 shown in Fig. 2.16.2, where n^, 
n 2 ,n^ are any mutually perpendicular unit vectors, one can formulate ^ as follows 
Let § 2 , a« be mutually perpendicular unit vectors directed as shown in 
Fig. 2.16.2. Define l^j^ as 

IjklSj*i‘Sk (j,k- 1,2,3) 
and form (i « 1,2,3) and V (R,X,3) as 


Qll A 9^v* 
3r2 


(9) 


( 10 ) 
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and 


Then 


Qoo 4 1 3V* . sec^S 9^V* tang 3V* 

2 2 “ R ^ ^ 

Q 33 A 1 3V* 1 32 ?* 


V (R,X,3) A mV - 2 (^1x^11 '*’ ^22^22 ^33^33 ) '*’ ^ 


F.o |I-+a ISM^+a JL|V 

- Si - 55 - + §2 R 3X ^ S 3 R jjf- 


( 11 ) 

( 12 ) 

(13) 

(14) 


For the derivation of Eq, (14), it is helpful to note that, since the 
indicated partial differentiation is referred to the reference frame in which 
n^, ^ 2,113 in Fig. 2.16.2 are fixed. 



3a^ 

2 Si 

w 


0 

c3a2 §3 

^§2 

3r” 

9®2 

W 

^S 2 

33 

ss 

0 

-c3aj^+s3a3 0 

9S3 

3R 

3a3 

3X 

3^3 

w 


0 

-s3a2 -a^ 


(15) 


Also, proceeding as in the Example in Sec . 2.9, one can verify that, if 


V(R) = V (R,X, 6 ), then 


3V . sec3 3V 1 3V 

S 2 R 3X S 3 R 33 


VV(R) = + §0 ^ + a, ^ 


(16) 


Differentiating once more with respect to R and using Eq. (15) one then 
finds that 


VVV(R) ““ Qxx«x^X ^22S2“2 ^33S3 — 3 


(17) 


where Q 22 » Q 33 given by Eqs. (4)- ( 6 ) while the dots represent all 

terms involving the dyads a^^a^ with i j. Similarly, the inertia dyadic I can 
be expressed as 
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^22-2^2 ^332323 + ••• 


Hence, 

I : ^11^11 ■*■ ^22^22 ■'■ ^33^33 

With V In place of V(R), one thus obtains from Eq. (2) 

V(R) . mV* - i + h3’>3l) * <= 

Moreover, defining V (R,X,@) as the right-hand member of this equations, which 
yields Eq. (13), one can write [compare with Eq. (16)] 


vm/n\ 3V . sec3 3V ^ 3V 

/ " 2i ^ + 22 ^ ^ + 23 ^ 

, / 

Subsfltu|tlon Into Eq. (1) then produces Eq. (14). 
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2,17 Force function expressions for moment exerted on a body by a particle 


Force functions (see Secs, 2,10 - 2.16) cap be used also in connection 
with gravitational moments. For the interaction of a particle and a (not 
necessarily rigid) body, the moment M given by Eq. (2.6.1) can be expressed as 
M = -R X VV (1) 

where V(R) is given by Eq. (2.11.2) or by Eq. (2,11.5) and VV is the deriva- 
tive of V with respect to R; or, if W is defined as W(a^^,R) ^ V[R(aj^,R)], M is 
given also by 


“ ”^i ^ '51^ 


( 2 ) 


Similarly, for M as defined in Eq, (2.6.3), one can write 

M = -R X VV (3) 

where V is given by Eq. (2.12,13); or M can be expressed as 

(4) 


M V 

M = -a- X ■!— 
1 9aj^ 


where W is a function of a^^ and R defined in terms of the moment of inertia 

of B about the line connecting P and B as 

. 3GmI.^ 

W = ^ + C 

2R'^ 


Derivations ; Eq. (1) is obtained by substituting from Eq. (2.11.1) into 
Eq. (2.6.1); and Eq. (2) then follows if one uses Eq. (2.9.12). 

Referring to Eqs. (2.3.6) and (2.3.3), one can write 
_ tr 3Gm 

-R a, X F =• — y a- X I . a, = M 

rJ -1 - -1 ( 2 . 6 . 3 ) ■ 


and, replacing R a^^ with R, one then arrives at Eq. (3) after using Eq. 
(2.12.12) to eliminate F. 


121 



2.17 


w 

If W is defined as 


W(aj^,R) A v[R(aj^,R)l 


then 


and 


W 


Gim& 


(2.12.13) 


1 + 


2 mR 

3W 


[tr(I)-3l3^J 


+ C 


R X W = a, X ^ 
" “ (2.9.12) 


- ii » H ^ * <=> 


or*, since tr (I) is independent of a^, 

3GmI, 


R X yv « a^^ X ^ 


2 l V 2r3 ) 


V 

” a - X 

(5) ®5l 


Consequently 


g . - a, X |» 

- (3) -I sjj 


in agreement with Eq. (4). 

Example , in Fig. 2.17,1, a^^, § 2 , and ^ 2 * b^ are dextral sets of’ 

orthogonal unit vectors, P is a particle of mass m, and B is a rigid body of 

mass m. The unit vector a^ is chosen such that the position vector R of the 
* _ 

mass center B of B relative to P is given by R = Ra^^; and bj^, b 2 * b^ are 
fixed in B. 


If 9j^, 02 , body-two orientation angles (see Sec. 1.7) for B in 

a reference frame in which a^^, § 2 , a^ are fixed, one can express R as 

R(R,ej^,02.02) " ^(®2-l ®2®3-2 ®2‘^3-3^ 

(1 . 7 . 31) 


Given a force function V(R) for the gravitational interaction of P and B, one 
can, therefore, define a function W(R, 0 j^, 02 » 63 ) as 


122 




2.17 


W(R, 0 ^, 02 , 63 ) 4 V[R(R,02,03)] 

* 


(7) 


The moment M exerted by P on B about B can then be expressed In terms of 
partial derivatives of W with respect to 03 and 03 by proceeding as follows: 
Partial differentiation of R with respect to R| 02 » and 03 yields 


3R 

3R *^2-1 ®2®3-2 ®2*^3-3 

(8) 

3R 

W 2 + ‘=2®3^2 +;‘^2‘^3^3> 

(9) 

3R 

^3 (6) *‘*2V3Si - '2®3S3> 

(10) 


From Eq. (2.9.5) 
3W 


3R 
^(8) 
3R 


tS - ^ • 5ii ■,‘=2 + » 2»3 ^ • ^2 + = 2=3 ^ • h 


(11) 


^ ‘ '^2 (9)* ^ • Si °2=3 2V • 62 + <=2'=3 2V • ^ 3 ) 


9W 

30, 


3R 


VV • 


W, ‘‘(= 2'=3 JV • bj - > 2=3 ^ • S 3 ) 


'3 " ”“3 ( 10 ) 

Eqs. (11) -(13) may be solved for the dot-products of VV with bj^, b 2 » ^ 3 ; and 
W can then be expressed as 

VV » VV • bj_b 3 + VV • 8382 + VV • 8383 


(13) 


-/Mo i 

y3R *^2 ~ R 302 2/ -1 

. /3W . 1 3W . 1 3W ®3\. 

\W ®2®3 R *'2®3 + R ^ "q j ^2 

. / 3W , 1 3W 1 3W ®3 \ , 

\3R ®2*^3 R 302 ®2'^3 " R 303 


(14) 
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Consequently 


M = -R X VV » 


(1) 


3W 

30, 


(6,14) “'^3 


h-( 


3W 

w" ‘=3 ■ m 


3W '^2®3 


3 ®2 




/3W 

C 2 C 3 

\90, 

So 


■*' 302 ®3/ -3 


}) h 


This equation does not Involve 0. explicitly; but, referred to a,, a, 

X —X — i 

chan to bj^,b 2 ,b 2 , M is given by {use (1,17.31)3. 

^ “ (‘^1 Hj) S2 + ^ s| if-) 23 

which brings the dependence of M on 0^ into evidence. 
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2.18 Force function expressions for moment exerted on a body by a body 

If V is the force function for a body B and a particle P of unit masSy 

* 

then the moment exerted by a body B on B about its mass center B is given by 

M = “ £ I ^ Pd^ (1) 

where,, as in Fig. 2.7.1, P is a generic point of 1, R is the vector from P 
to B , p is the mass density of B and dx is a differential element of volume 
of B. 

9tm 

When the distance R between the mass centers B and B of B and B is 
large relative to the largest dimensions of B, then M can be approximated by 


pdT 

B ^ 

in which R and a^^ are defined by R = R a^^, as in Fig. 2.7.1'. 

When R is large relative to the dimensions of both B and B, M can be 
approximated by Eqs. (2.17.3) or Eqs. (2.17.4) and (2.17.5) with V from Eqs. 
(2.12.13), where now R, R, and a^^ are to be interpreted as indicated in Fig. 

"k 

2.7.1, is the moment of inertia of B about the line joining B and B , 
and m is the mass of B. 

Derivations t 

"it 

Eq. (1) follows from the definition of a moment about B and the identi*- 
fi cation of VgV pdx as the force exerted on B by the differential element of 
mass at P. 

As shown in Eqs. (2.12.1) and (2.12.2), the force function for B and an 
element of unit mass at P can be written as 


V 


Gm ( , , 1 

* ~ ^ + ~T2 

R ( 2mR 






(3) 
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where represents terms of degree in the ratio jrl/R (see Fig* 2.7.1). 

Differentiation yields 



Substitution of Eq. (4) into Eq. (1) leads to 



R = R(1 + ...) 

(6) 

a^^ * a^^ . 

(7) 


where three dots represent terms of first degree or higher in the ratio 
|rl/R (see Fig. 2.7.1). By substituting Eqs. (6) and (7) into Eq. (2), and 
using the definition 

m — ^ pdx (8) 

one can confirm the applicability of Eq. (2.7.2), with a^^ now interpreted 
as in Fig. 2.7,1; Eq. (2,7.2) has been shown in Section 2.17 to be equivalent 
to Eqs. (2*17.3) and to the combination of Eqs. (2.17.4) and (2.17,5). 
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Example 

Fig. 2.4.2 In Section 2.4 portrays a homogeneous oblate spheroid B and 
a homogeneous rectangular parallelepiped with the mass centers B and B of 
B and B respectively separated by a distance R that substantially exceeds the 
maximum dimensions of B and B« 

mm . ife 

The moment exerted by B on B about B Is approximated by 

M . - R X V V (9) 

(2.17.3) " 


where 


, o ji + -L_ ftrd) - 3a, • I • a,l | 

(2.12.13) ^ ( 2mR^ L ^ ■ J ) 

As Implied by the Example of Section 2.14, Eq. (10) produces 

2 2 


( 10 ) 




+ G 


( 11 ) 


[This result should be contrasted with V as it actually appears In the Example 
of Section 2.14; the same distinction exists between the V of Eq. (2.12.13) and 
the V of Eq. (2.14.3)]. 

Differentiation of Eq. (11) provides 


V_V 

- (11,2.9.8, 
2.9.9) 


Gmm 

V 


2 2 r 

*1 i 0 \ 

^ CC e 

j|^l-3Cii) a^ + 


( 12 ) 


as indicated by the Examples in Sections 2.14 and 2.4. 

Substitution and cross multiplication produce 

*-2 2 / 

M » Gmma e ( C 

"( 9 , 12 ) ^^3 ^ 

This result should be Identical to Eq. (2.7.6) in the Example of 
Section 2.7. To reconcile these two expressions for M, one must recognize 
the identities 


21-3 


- C 


3i: 


22) 


(13) 
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I- - * ma^e^/12 

/(2.7.9, 

2.7.10, 

2.4.25) 


and 

^3-2 “ ^12-3 “ ^21-3 ~ ^31®2 
The validity of Eq. (15) follows from the expressions 

-3 “ ^31-1 ^32-2 ^33-3 

-2 ” ^21-1 *'22-2 *'23-3 


(14) 


(15) 


(16) 

(17) 


and the equality of an element of the direction cosine matrix with Its 
cofactor, both of which are established In Section 1.2. 
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2.19 Force function expression for the moment exerted on a small body by a 
body 

When the distance between the mass centers B and B of two bodies B 
and B exceeds the greatest distance from B to any point of B, and a force 
function V(g) for the forces exerted by B on a particle of unit mass at a 
point F situated as shown In Fig. 2.16.1 Is available, then the system of 

* A 

gravitational forces exerted by B on B produces a moment about B that is 
given approximately by 

M - -I X V7V(R) (1) 

where I Is the Inertia dyadic of B for B , R Is the position vector of B 

-* t 

relative to B , and V connotes differentiation with respect to R. 

Derivation; The system of gravitational forces exerted by B on B 


< 2 ) 


produces a moment M about B such that 

M = J't ^ 

Expanding V^V in a Taylor series as in Sec. 2.16 following Eq. (2.16.4), one 
can sho\^ that 

V^V = VV(R) -f r • WV(R) ... 

so that, retaining only terms displayed explicitly, one arrives at 

M * Tr X VV(R) pdx + y*r X |^r • VW(R)1 pdT 

( 2 ) 

* 

The first integral vanishes because r originates at the mass center of B ; 
and 

r X ^r • WV(R)j “ r r x VW(R) 


The cross-dot product in Eq. (1) is defined such that, for two dyads 
u^u^ - 1 Z 2 * • ^-1~2^ "" ^-2* -2^* and it obeys the dis- 

tributive law when applied to dyadics: 

(a2^a2"f"bj^b2"H . • « ) x . • . ) — ^1^2^^1^2"^^1^2^^1^2'^ • • • 

«1-2*~1-2^-1-2«-i52^* * * 

+ • . . 
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Hence 

j tr(I)U-lj X WV(R) 

But WV(R) is a symmetric dyadic^ and the cross-dot product of U and any 
symmetric dyadic is equal to zero® What remains is Eq® (1)® 

Example ; When V(R) is available in the form of an explicit function V (R,X,3) 
of the spherical coordinates R,X>3 shown in Fig. 2.16.2, one can find M as 
follows : 

Differentiate Eq. (2.16.16) to verify that 

WV(R) = ^ i Q±a a. a 
i==l j=i ^ J 


IX 

M 




pdT ] X WV(R) = 


where Q22* ^33 given in Eqs. (2.16.10 - 2.16.12) and 



Next , let 

= fj * J * 2k (j,k=l,2,3) 

Then 

“ ” *^12^31 ^^22"^33^^23 " ^^22~^33^^23]-l 

[^23^12 " ^23^12 ^^33"^11^^31 " ^^33"^11^^3l]-2 

[^31^23 ~ ^31^23 ^^ll"^22^^12 “ ^^ll"^22^^12]-3 
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2.20 Gravitational potential energy 

In the context of Lagranglan ot Hspiltonlan mechanics, the phrase 
"potential energy” refers to a function of the generalized coordinates 
(r-l,...,n) of a system S (in an inertial reference frame A) such that 


3* 

^ (r*l,...,n) 


( 1 ) 


where F^, the generalized force corresponding to q^, is defined in terms of 

the position vector of a typical particle of S relative to a point fixed 

in A, and F^, the resultant of all contact and body forces acting on as J 

A ^ 3P. ' ' 

A V* -i 

i«l *^ ^r 


F A £ ^ . 

r fr, 3q,, -i 


(r»l, , . , ,n) 


m 


Here N is the number of particles and the partial differentiations are perfoimed 
in reference frame A. In some situations, $ is simply the negative of one of 
the force functions discussed in Secs. 2.10 - 2.19. For example, this is tito 
case when S consists of a single particle P of mass m moving under the action 
of the gravitational force exerted by a particle P of mass m fixed in A. $ is 
then given by 


$ ss _ (Gmmp ^ + C) (3) 

where G, p, and C have precisely the same meaning as in Eq. (2.10.2). Eq, 

(3) also furnishes the potential energy when S consists of two particles mov'- 
Ing under the action of their mutual gravitational attraction. In other 
situations, $ cannot be found so simply, but some of the force functions of 
Secs. 2.10-2.14 are, nevertheless, useful, for they provide contributions • 
to $. For example, if S consists of a rigid body B of mass m and two particles 
P^ and p 2 , of masses m^ and m 2 , respectively, and if these three objects move 


The symbol # is used for potential energy rather than the more familiar V, 
to avoid confusion with the use of V for the force function in Secs. 2.10-2.19. 
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under the action of forces of mutual gravitational attraction, then $ can be 
expressed as 

$ = 

with [see Eq. (2.10.2)] 

$ » Gm_m„p”^ + C (5) 

o 1 2*^0 o 

and [see Eq. (2.11.5)] 

^1 “ * ^1 (l-lrZ) (6) 

where p Is the distance between P, and ?o, Pj Is the distance from P. to 
a generic point P of B, p Is the mass density of B at P, and C^, and C 2 are 

arbitrary constants. 

Force functions associated with approximations to gravitational forces 
and moments can at times be used In the formulation of a potential energy 
but ,lt can occur that apparently relevant force functions exist when $ does 
not exist. For example, let S consist of a rigid body B of mass m moving under 
the action of gravitational forces exerted by a particle P of mass m fixed In 
A, and suppose that the resultant gravitational force F exerted on B by P Is 
approximated with F as defined In Eq. (2.3.6), while the moment M about the 

is <%i# 

mass center B of B of all forces exerted on B by F Is approximated with M as 
defined In Eq. (2.6.3). Then the potential energy $ of S can be expressed as 
$ . -V (7) 

with V given by Eq, (2,12#13); but if ? is approximated with F as defined in 

rw 

Eq« (2e3*5)^ while M is approximated with then (as shown for a specific 

case in the Example following) there exists no function ^ that satisfies Eq. 

'A 

(1), despite the fact that force functions for F and H do exist in the sense 
A rs# 

that F and M can be expressed as 
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F - («> 

and 

M - -R X (9) 

ft M 

where R Is the position vector of B relative to F while 

- GmnR"^ (10) 

and 

V 2 ■ G5mR”^v^^^ (11) 

( 2 ) 

with V ' as given in Eq. (2.12.2). (The non-existence of a potential energy 
under these circumstances may be an Inconvenience, but it does not render the 
approximations under consideration invalid. They are, in fact, particularly 
useful, and are, therefore, employed extensively, in the analysis of motions 
of space vehicles). 

Derivations ; The resultant of all contact and body forces acting on a 
single particle P of mass m moving under the action of the gravitational force 
exerted by a particle P of mass ffl is given by 

F = -GSnp(p^)”^^^ (12) 

" ( 2 . 1 . 1 ) 


Consequently 


^r " 
( 2 ) 



3E 2 ' 

F ■ -Gmm ■ 5 — ♦ p (p ) 

■ (12) ®‘*r - - 



■ T— 6 mmp‘ 
3qj. 



which shows that $ as given in Eq. (3) is the potential energy of a system S 
consisting solely of P. Similarly, if S consists of both P and P, the 
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When S consists of a rigid body B of mass m and two particles, and 
P 2 * masses m^ and m 2 , the resultants and £2 of all contact and body forces 
acting on and ?2 are respectively [see Eqs. (2.1.1), (2.1.2) ,(2.2.2)] 

Pj - -Cjmj 2o(£)'^'^ * ®"i/ Ei(Ei) 

and 

h ■ + ®"2 Jt2 ( 

where p^ Is the position vector of m^^ relative to m 2 , and p^ and P 2 are the 
position vectors of a generic point P of B relative to P^ and ?2 as shown In 
Fig. 2.20.1. The contribution to of forces acting on P^ and P 2 is given 
by [see Fig. 2.20.1 for and R 2 , and use Eq. (2)]. 
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(17) 


Ai for forces acting on the particles of B, these either are forces exerted 
by iui^ P2> or they are forces exerted by particles of B on each other. 

$ince B Is a rlBld body, the latter forces make no contribution to (see Kane, 
bvnaalcs. Holt, Rinehart and Winston, p. 81), and a typical force of the former 
kind can be expressed as 


r ! 

2V-3/2 , 



dF - -6 m-p-j p 

(2.1.1) L 

l) "2?2(? 

2 ) 

(18) 


The contribution of this force to Is thus given by [see Fig, 2.20.1 for g, 
and use Eq. (2)] 



And can now be formulated as 
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9R 9R ^ 3p 

• F, + -5-^ « F, + I -5^ 


(17,19) ^ ^ ^**1 


?o(Eo) 


- /■ ^^2 / 2 
y ^ • ?l(?l) l='‘^ - “27 3|; • f 2 (E 2 ) <>'*■' 


[®"i'“2(Eo) ^ +°”i/(bJ) 


pdT + Gm. 




A / 2\^^^ 

“(^i) 

= 9^1' 

(20) L 


(j = 0 , 1 , 2 ) 


Gm^itt2 pQ^ ■*■ 0®2^ ^ pdx + 0®2^ ^2^ 


(5,6) 


(*0-^*1 -^* 2 ) 


(4) 


Hence, $ as given by Eq. (4) is the potential energy of S. 

To form the generalized force associated with the force system presumed 
to be acting on a rigid body B in connection with Eq. (17) , one can begin by 

expressing F^ as (see Kane , . Dynamics , Holt, Rinehart and Winston, p. 81) 

~ a , > 9o) 

where R and a^ are defined as in Sec. 2.13, while w, the angular velocity of 
B in reference frame A, is regarded as a function of and (r-1, . . . ,.n), . 
Using Eqs. (2.3.6) and (2.6.3), one thus obtains 


^ — Gmm 

F = 

^(21) R'^ 




or, equivalently. 
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(22) R 

4 . 

r3 


Gmm 

2 1-1 3q 


* K) (■^i) 


" Si)- (i • 2i) 


Now, 


3a^ 


3 / „ \ 3R . _ ^ 

3q^ (®2l) “ 34^ 2i + R 3q 


and, using a superscript B to denote differentiation In B, one can write 

3i;“ ^ + 


so that 


9 a, 


j(25) ^9 

- • I 
‘r " 

m 

a, - ^ • I 

-1 3qj. - 

q (r=l 
^r 

J • • • > 

n) 

In B, 

\ 

9 Si 


1®3 

-l) “ 

^'(26) 

w * 

I 

* -1 ■ 2 3q^ 


^—1 

T 

1 3^11 



X 

-l" 2 3q 

r 


where is the moment of inertia of B about the line passing through B 
and parallel to a^^. Substituting from Eqs. (24) and (27) into Eq, (23), 
thus arrives at 


„ + Ba . ^ + ^ a • f ^2) + R !|l . f (2)1 

®^r“ ^2|^3q^ + »Si 3q^ + 3q^ Si I + » 3q^ f J 

. 3^ / 1 3^ii\ 

\3 i^3q/l*2l-2 3qj 

The vector 3aj^/3q^ Is necessarily perpendicular to a^^. Consequently the 

( 2 ) 

dot-product of this vector with a^^ vanishes. Expressing f as In Eq, 


(23) 

(24) 

(25) 

(26) 

(27) 

I; 

one 

(28) 
(28) 

(2.3.3), 
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one can thus rewrite Eq. (28) as 

1 3R L Ij. 9^11 

ltr(I)-3I,J+ ^ 

3V 

" W 
(2.12.2) ^ 

which shows that $ as given in Eq. (7) is the potential energy in question. 

Example; Fig. 2,20.2 shows a thin rod B of mass m and a particle P of 
mass m, the same system previously cbnsldered in the Examples in Sec, 2.3 and 
2.11. If B is constrained to remain in a plane fixed in an inertial refer- 
ence frame A, and if P is fixed in this plane, then the position and orientation 
of B in A can be described in terms of the generalized coordinates q^ 6, 
q 2 = R, and q^ = ip. Suppose now that F, the resultant force exerted by P on B, 
is approximated as 


Gmm 3R 3Gm 

2 ” 3 

R 3q^ 2 r'* 


while M, 
with 


FssF - -^^fa- + -^ (2-3 sin^^p) a. - 

" ■(2.3.6) ^ L"" 8r2 ^ ' 8r2 

it . 

the moment about B of the forces exerted on B by P, 


M ssM - - ^^(sin2\p a-) 

" (2,6.3) 8R-^ ' 


sin2ip a 2 j (29) 
is approximated 

(30) 


Then the generalized forces corresponding to 6,R,ip are 



(31) 

(32) 

(33) 
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Figure 2.20.2 
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and, if $ is taken to be the negative of the function V found in the Example 
in Sec. 2.11, that is. 


$ = > 


Gmm I 1 + 


(2-3 sin^i|>) 


R 1 24R 

Then, as may be verified by carrying out the indicated differentiations. 


t? _ ^ 1? _ ^ If 

" " 30 » ^R ■ 3R’ “ " 3i|> 

But, if instead of approximating F as in Eq. (29) one uses 

_ A gum 

F « F = - -r~ a, 

R^ -J. 

while continuing to approximate M as in Eq. (30) , then the generalized 
forces become 


(34) 


(35) 


'e ■ 'n, ■ - 


GnnnL ^ . « , « Gmm 

.1.2^ F, - - 


and there exists no function $ of 6, R and if) such that Eqs. (34) are satisfied 

simultaneously. Finally, if additional constraints are Introduced so as to 

force B to move with a prescribed velocity, which can be accomplished by 

A 

applying suitable contact forces at B , then 9 and R cease to be generalized 
coordinates! the generalized force for what is now a one-degree-of -freedom 
system, is given by Eq. (33), both when F is approximated as in Eq. (29) and 
when Eq. (35) is used; and a potential energy $ thus exists in both cases 
because F^ can be expressed as F^ = - if one lets 

» = . + c 

16R 
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